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Description 

BIOINFORMATICALLY DETECTABLE 
GROUP OF NOVEL REGULATORY 
OLIGONUCLEOTIDES ASSOCIATED WITH 
ALZHEIMER'S DISEASE AND USES 

THEREOF 

Cross Reference to Related Applications 

[0001] This application is a continuation in part of and claims 

priority from the following patent applications, the disclo- 
sures of which applications are all hereby incorporated 
herein by reference: U.S. Patent Application Serial No. 
10/707975 filed 29-Jan-04,U.S. Patent Application Serial 
No. 10/707147 filed 24-Nov-03, U.S. Patent Application 
Serial No. 10/707147 filed 24-Nov-03, U.S. Patent Appli- 
cation Serial No. 10/604985 filed 29-Aug-03, U.S. Patent 
Application Serial No. 10/651227 filed 29-Aug-03, U.S. 
Patent Application Serial No. 10/649653 filed 28-Aug-03, 
U.S. Patent Application Serial No. 10/604926 filed 



27-Aug-03, U.S. Patent Application Serial No. 10/604726 
filed 13-Aug-03, U.S. Patent Application Serial No. 
10/604727 filed 13-Aug-03, and U.S. Provisional Patent 
Application Serial No. 60/468251 filed 07-May-03. This 
application also claims priority from International Applica- 
tion Number: PCT/IL 03/00970, filed 16-Nov-03, the dis- 
closure of which application is hereby incorporated herein 
by reference. All of the aforesaid patent applications are 
entitled "Bioinformatically Detectable Group of Novel Reg- 
ulatory Genes and Uses Thereof"; This application also 
claims priority from U.S. Patent Application Serial No. 
10/707980 filed 29-Jan-04,entitled "Bioinformatically De- 
tectable Group of Novel Regulatory Oligonucleotides and 
Uses Thereof; U.S Patent Application Serial No. 
10/707980, filed 29-Jan-04, entitled "Bioinformatically 
Detectable Group of Novel Regulatory Oligonucleotides 
and Uses Thereof" is a continuation in part of and claims 
priority from the following patent applications, the disclo- 
sures of which applications are all hereby incorporated 
herein by reference: U.S. Patent Application Serial No. 
10/707147 filed 24-Nov-03, U.S. Patent Application Serial 
No. 10/707147 filed 24-Nov-03, U.S. Patent Application 
Serial No. 10/604985 filed 29-Aug-03, U.S. Patent Appli- 



cation Serial No. 10/651227 filed 29-Aug-03, U.S. Patent 
Application Serial No. 10/649653 filed 28-Aug-03, U.S. 
Patent Application Serial No. 10/604926 filed 27-Aug-03, 
U.S. Patent Application Serial No. 10/604726 filed 
13-Aug-03, U.S. Patent Application Serial No. 10/604727 
filed 13-Aug-03, and U.S. Provisional Patent Application 
Serial No. 60/468251 filed 07-May-03. This application 
also claims priority from International Application Num- 
ber: PCT/IL 03/00970, filed 16-Nov-03, the disclosure of 
which application is hereby incorporated herein by refer- 
ence. All of the aforesaid patent applications are entitled 
"Bioinformatically Detectable Group of Novel Regulatory 
Genes and Uses Thereof; U.S Patent Application Serial No. 
10/707975, filed 29-Jan-04, entitled "Bioinformatically 
Detectable Group of Novel Regulatory Genes and Uses 
Thereof" is a continuation in part of and claims priority 
from the following patent applications, the disclosures of 
which applications are all hereby incorporated herein by 
reference: U.S. Patent Application Serial No. 10/707147 
filed 24-NOV-03, U.S. Patent Application Serial No. 
10/707147 filed 24-Nov-03, U.S. Patent Application Serial 
No. 10/604985 filed 29-Aug-03, U.S. Patent Application 
Serial No. 10/651227 filed 29-Aug-03, U.S. Patent Appli- 



cation Serial No. 10/649653 filed 28-Aug-03, U.S. Patent 
Application Serial No. 10/604926 filed 27-Aug-03, U.S. 
Patent Application Serial No. 10/604726 filed 13-Aug-03, 
U.S. Patent Application Serial No. 10/604727 filed 
13-Aug-03, and U.S. Provisional Patent Application Serial 
No. 60/468251 filed 07-May-03. This application also 
claims priority from International Application Number: 
PCT/IL 03/00970, filed 16-Nov-03, the disclosure of 
which application is hereby incorporated herein by refer- 
ence. All of the aforesaid patent applications are entitled 
"Bioinformatically Detectable Croup of Novel Regulatory 
Genes and Uses Thereof; U.S Patent Application Serial No. 
10/707147, filed 24-Nov-03, entitled "Bioinformatically 
Detectable Group of Novel Regulatory Genes and Uses 
Thereof is a continuation in part of and claims priority 
from the following patent applications, the disclosures of 
which applications are all hereby incorporated herein by 
reference: U.S. Patent Application Serial No. 10/604985 
filed 29-Aug-03, U.S. Patent Application Serial No. 
10/651227 filed 29-Aug-03, U.S. Patent Application Se- 
rial No. 10/649653 filed 28-Aug-03, U.S. Patent Applica- 
tion Serial No. 10/604926 filed 27-Aug-03, U.S. Patent 
Application Serial No. 10/604726 filed 13-Aug-03, U.S. 



Patent Application Serial No. 10/604727 filed 13-Aug-03, 
and U.S. Provisional Patent Application Serial No. 
60/468251 filed 07-May-03. This application also claims 
priority from International Application Number: PCT/IL 
03/00970, filed 16-Nov-03, the disclosure of which ap- 
plication is hereby incorporated herein by reference. All of 
the aforesaid patent applications are entitled "Bioinformat- 
ically Detectable Group of Novel Regulatory Genes and 
Uses Thereof; International Application Number: PCT/IL 
03/00970, filed 16-Nov-03, entitled "Bioinformatically 
Detectable Group of Novel Regulatory Genes and Uses 
Thereof" is a continuation in part of and claims priority 
from the following patent applications, the disclosures of 
which applications are all hereby incorporated herein by 
reference: U.S. Patent Application Serial No. 10/604985 
filed 29-Aug-03, U.S. Patent Application Serial No. 
10/651227 filed 29-Aug-03, U.S. Patent Application Se- 
rial No. 10/649653 filed 28-Aug-03, U.S. Patent Applica- 
tion Serial No. 10/604926 filed 27-Aug-03, U.S. Patent 
Application Serial No. 10/604726 filed 13-Aug-03, U.S. 
Patent Application Serial No. 10/604727 filed 13-Aug-03, 
U.S. Provisional Patent Application Serial No. 60/468251 
filed 07-May-03, and U.S. Patent Application Serial No. 



10/345201 filed 16-Jan-03. All of the aforesaid patent 
applications are entitled "Bioinformatically Detectable 
Group of Novel Regulatory Genes and Uses Thereof; U.S 
Patent Application Serial No.10/604985, filed 29-Aug-03, 
entitled "Bioinformatically is a continuation of U.S Provi- 
sional Patent Application Serial No. 60/468251, filed 
07-May-03, entitled "Bioinformatically Detectable Group 
of Novel Regulatory Genes and Uses Thereof " the disclo- 
sure of which is hereby incorporated herein and claims 
priority therefrom; and is a continuation in part of and 
claims priority from the following patent applications, the 
disclosures of which applications are all hereby incorpo- 
rated herein by reference: U.S. Patent Application Serial 
No. 10/651227 filed 29-Aug-03, U.S. Patent Application 
Serial No. 10/649653 filed 28-Aug-03, U.S. Patent Appli- 
cation Serial No. 10/604926 filed 27-Aug-03, U.S. Patent 
Application Serial No. 10/604726 filed 13-Aug-03, U.S. 
Patent Application Serial No. 10/604727 filed 13-Aug-03, 
U.S. Patent Application Serial No. 10/345201 filed 
16-Jan-03, U.S. Patent Application Serial No. 10/321503 
filed 18-Dec-02, U.S. Patent Application Serial No. 
10/310914 filed 06-Dec-02, and U.S. Patent Application 
Serial No. 10/293338 filed 14-Nov-02. All of the afore- 



said patent applications are entitled "Bioinformatically De- 
tectable Group of Novel Regulatory Genes and Uses 
Thereof"; U.S Patent Application Serial No. 10/604926, 
filed 27-Aug-03, entitled "Bioinformatically Detectable 
Group of Novel Regulatory Genes and Uses Thereof is a 
continuation of U.S Patent Application Serial No. 
10/345201, filed 16-Jan-03, entitled "Bioinformatically 
Detectable Group of Novel Regulatory Genes and Uses 
Thereof" the disclosure of which is hereby incorporated 
herein and claims priority therefrom; and is a continuation 
in part of and claims priority from the following patent 
applications, the disclosures of which applications are all 
hereby incorporated herein by reference: U.S. Patent Ap- 
plication Serial No. 10/604726 filed 13-Aug-03, U.S. 
Patent Application Serial No. 10/604727 filed 13-Aug-03, 
U.S. Provisional Patent Application Serial No. 60/468251 
filed 07-May-03, U.S. Patent Application Serial No. 
10/321503 filed 18-Dec-02, U.S. Patent Application Serial 
No. 10/310914 filed 06-Dec-02, and U.S. Patent Applica- 
tion Serial No. 10/293338 filed 14-Nov-02. All of the 
aforesaid patent applications are entitled "Bioinformati- 
cally Detectable Group of Novel Regulatory Genes and 
Uses Thereof; U.S Patent Application Serial No. 



10/649653, filed 28-Aug-03, entitled "Bioinformatically 
Detectable Group of Novel Regulatory Genes and Uses 
Thereof" is a continuation of U.S Patent Application Serial 
No. 10/321503, filed 18-Dec-02, entitled "Bioinformati- 
cally Detectable Group of Novel Regulatory Genes and 
Uses Thereof "; the disclosure of which is hereby incorpo- 
rated herein and claims priority therefrom; and is a con- 
tinuation in part of and claims priority from the following 
patent applications, the disclosures of which applications 
are all hereby incorporated herein by reference: U.S. 
Patent Application Serial No. 10/604926 filed 27-Aug-03, 
U.S. Patent Application Serial No. 10/604726 filed 
13-Aug-03, U.S. Patent Application Serial No. 10/604727 
filed 13-Aug-03, U.S. Provisional Patent Application Serial 
No. 60/468251 filed 07-May-03, U.S. Patent Application 
Serial No. 10/321503 filed 18-Dec-02, U.S. Patent Appli- 
cation Serial No. 10/310914 filed 06-Dec-02, and U.S. 
Patent Application Serial No. 10/293338 filed 14-Nov-02. 
All of the aforesaid patent applications are entitled "Bioin- 
formatically Detectable Group of Novel Regulatory Genes 
and Uses Thereof"; U.S Patent Application Serial 
No.10/651227, filed 29-Aug-03, entitled "Bioinformati- 
cally Detectable Group of Novel Regulatory Genes and 



Uses Thereof is a continuation of U.S Patent Application 
Serial No. 10/310914, filed 06-Dec-02, entitled "Bioinfor- 
matically Detectable Group of Novel Regulatory Genes and 
Uses Thereof"; the disclosure of which is hereby incorpo- 
rated herein and claims priority therefrom; and is a con- 
tinuation in part of and claims priority from the following 
patent applications, the disclosures of which applications 
are all hereby incorporated herein by reference: U.S. 
Patent Application Serial No. 10/604985 filed 29-Aug-03, 
U.S. Patent Application Serial No. 10/649653 filed 
28-Aug-03, U.S. Patent Application Serial No. 10/604926 
filed 27-Aug-03, U.S. Patent Application Serial No. 
10/604726 filed 13-Aug-03, U.S. Patent Application Se- 
rial No. 10/604727 filed 13-Aug-03, U.S. Provisional 
Patent Application Serial No. 60/468251 filed 07-May-03, 
U.S. Patent Application Serial No. 10/345201 filed 
16-Jan-03, U.S. Patent Application Serial No. 10/321503 
filed 18-Dec-02, U.S. Patent Application Serial No. 
10/310914 filed 06-Dec-02, and U.S. Patent Application 
Serial No. 10/293338 filed 14-Nov-02. All of the afore- 
said patent applications are entitled "Bioinformatically De- 
tectable Group of Novel Regulatory Genes and Uses 
Thereof"; U.S Patent Applications Serial Nos. 10/604727 



and 10/604726, filed 13-Aug-03, entitled "Bioinformati- 
cally Detectable Group of Novel Regulatory Genes and 
Uses Thereof are a continuation of U.S Patent Application 
Serial No.10/293338, filed 14-Nov-02, entitled "Bioinfor- 
matically Detectable Group of Novel Regulatory Genes and 
Uses Thereof the disclosure of which is hereby incorpo- 
rated herein and claims priority therefrom; and are a con- 
tinuation in part of and claims priority from the following 
patent applications, the disclosures of which applications 
are all hereby incorporated herein by reference: U.S. Pro- 
visional Patent Application Serial No. 60/468251 filed 
07-May-03, U.S. Patent Application Serial No. 10/345201 
filed 16-Jan-03, U.S. Patent Application Serial No. 
10/321503 filed 18-Dec-02, U.S. Patent Application Serial 
No. 10/310914 filed 06-Dec-02, and U.S. Patent Applica- 
tion Serial No. 10/293338 filed 14-Nov-02. All of the 
aforesaid patent applications are entitled "Bioinformati- 
cally Detectable Group of Novel Regulatory Genes and 
Uses Thereof; U.S Provisional Patent Application Serial 
No. 60/468251, filed 07-May-03, entitled "Bioinformati- 
cally Detectable Group of Novel Regulatory Genes and 
Uses Thereof is a continuation in part of and claims pri- 
ority from the following patent applications, the discio- 



sures of which applications are all hereby incorporated 
herein by reference: U.S. Patent Application Serial No. 
10/345201 filed 16-Jan-03, U.S. Patent Application Serial 
No. 10/321503 filed 18-Dec-02, U.S. Patent Application 
Serial No. 10/310914 filed 06-Dec-02, and U.S. Patent 
Application Serial No. 10/293338 filed 14-Nov-02. All of 
the aforesaid patent applications are entitled "Bioinformat- 
ically Detectable Group of Novel Regulatory Genes and 
Uses Thereof; U.S Patent Application Serial No. 
10/345201, filed 16-Jan-03, entitled "Bioinformatically 
Detectable Group of Novel Regulatory Genes and Uses 
Thereof" is a continuation in part of and claims priority 
from the following patent applications, the disclosures of 
which applications are all hereby incorporated herein by 
reference: U.S. Patent Application Serial No. 10/321503 
filed 18-Dec-02, U.S. Patent Application Serial No. 
10/310914 filed 06-Dec-02, and U.S. Patent Application 
Serial No. 10/293338 filed 14-Nov-02. All of the afore- 
said patent applications are entitled "Bioinformatically De- 
tectable Group of Novel Regulatory Genes and Uses 
Thereof"; U.S Patent Application Serial No. 10/321503, 
filed 18-Dec-02, entitled "Bioinformatically Detectable 
Group of Novel Regulatory Genes and Uses Thereof is a 



continuation in part of and claims priority from tlie fol- 
lowing patent applications, the disclosures of which appli- 
cations are all hereby incorporated herein by reference: 
U.S. Patent Application Serial No. 10/310914 filed 
06-Dec-02, and U.S. Patent Application Serial No. 
10/293338 filed 14-Nov-02. All of the aforesaid patent 
applications are entitled "Bioinformatically Detectable 
Croup of Novel Regulatory Genes and Uses Thereof; U.S 
Patent Application Serial No. 10/310914, filed 
06-Dec-02, entitled "Bioinformatically Detectable Group 
of Novel Regulatory Genes and Uses Thereof is a continu- 
ation in part of U.S Patent Application Serial 
NolO/293338, filed 14-Nov-02, entitled "Bioinformati- 
cally Detectable Group of Novel Regulatory Genes and 
Uses Thereof", the disclosure of which is hereby incorpo- 
rated by reference and claims priority therefrom. 
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Background of Invention 



FIELD OF THE INVENTION 



[0025] The present invention relates to a group of bioinformati- 

cally detectable novel oligonucleotides, here identified as 

Genomic Address Messenger or GAM oligonucleotides, 

which are believed to be related to the micro RNA (miRNA) 

group of oligonucleotides. 
DESCRIPTION OF PRIOR ART 

[0026] Micro RNAs (miRNA), are short ~22nt non-coding regula- 
tory RNA oligonucleotides, found in a wide range of 
species, believed to function as specific gene translation 
repressors, sometimes involved in cell-differentiation. 

[0027] The ability to detect novel mlRNAs is limited by the 

methodologies used to detect such oligonucleotides. All 
miRNAs identified so far either present a visibly discern- 
able whole body phenotype, as do Lin-4 and Let-7 
(Wightman,B., Ha,!., and Ruvkun,G., Cell 75:855-862 
(1993); Reinhart et al. Nature 403: 901-906 (2000)), or 
produce sufficient quantities of RNA so as to be detected 
by the standard molecular biological techniques. 

[0028] Studies reporting miRNAs (Lau et al.. Science 

294:858-862 (2001), Lagos-Quintana et al.. Science 294: 
853-858 (2001)) discovered 93 miRNAs in several 



species, by sequencing a limited number of clones (300 
by Lau and 100 by Lagos-Quintana) of small segments 
(i.e. size fractionated) RNA. MiRNAs detected in these 
studies therefore, represent the more prevalent among 
the mlRNA oligonucleotide family, and can not be much 
rarer than 1% of all small ~20nt-long RNA oligonu- 
cleotides. 

[0029] jhe aforesaid studies provide no basis for detection of 

miRNA oligonucleotides which either do not present a vi- 
sually discernable whole body phenotype, or are rare (e.g. 
rarer than 0.1% of all size fractionated ~20nt-long RNA 
segments expressed in the tissues examined), and there- 
fore do not produce significant enough quantities of RNA 
so as to be detected by standard biological techniques. 

[0030] Previous studies on mlRNAs and their relation to diseases 
have suggested potential involvement of several miRNAs 
in various type of cancers; It has been suggested that mir- 
15 and mir-16 are associated with B-cell chronic lympho- 
cytic leukemia (Calin,G.A at al., Proc. Natl. Acad. Sci. U. S. 
A, 2002). More recently, researchers have shown strong 
evidence for involvement of mir-143 and mir-145 in col- 
orectal neoplasia (Michael, M.Z. et al., Mol. Cancer Res. 1: 
882-891 (2003)). Mietzler and colleagues have demon- 



strated that mir-155, which is located on BIC locus, is 
highly and differentially expressed in pediatric Burkit lym- 
phoma patients (Metzler,M. at al. . Cancer 39: 167-169 
(2004)). Involevemnet of mlRNAs in Alzheimers disease is 
unknown. 

[0031] The following U.S. Patents relate to bioinformatic detec- 
tion of genes: U.S Patent No. 6,369,195, entitled "Prostate- 
specific gene for diagnosis, prognosis and management of 
prostate cancer", and U.S Patent No. 6, 291, 666 entitled 
"Spike tissue-specific promoter", each of which is hereby 

incorporated by reference herein. 
BRIEF DESCRIPTION OF SEQUENCE LISTING, LARGE TABLES 
AND COMPUTER PROGRAM LISTING 

[0032] A sequence listing is attached to the present invention, 
comprising 7,351 genomic sequences, is contained in a 
file named SEQ_LIST.TXT (1,207KB, 16-Feb-04), and is 
hereby incorporated by reference herein. 

[0033] Large tables relating to genomic sequences are attached 
to the present application, appear in 11 table files (size, 
creation date), incorporated herein: TABLE1.TXT (215 KB, 
16-Feb-04); TABLE2.TXT (16,358 KB, 16-Feb-04); TA- 
BLE3.TXT (134 KB, 16-Feb-04); TABLE4.TXT (945 KB, 
16-Feb-04), TABLE5.TXT (113 KB, 16-Feb-04), TA- 



BLE6.TXT (189 KB, 16-Feb-04) and TABLE7.TXT (3,335 KB, 
16-Feb-04), TABLE8.TXT (12,240 KB, 16-Feb-04), TA- 
BLE9.TXT (34,018 KB, 16-Feb-04), TABLE10.TXT (1,300 
KB, 16-Feb-04) and TABLEll.TXF (2 KB, 16-Feb-04), all 
of which are incorporated by reference herein. 
[0034] A computer program listing of a computer program con- 
structed and operative in accordance with a preferred em- 
bodiment of the present invention is enclosed on an elec- 
tronic medium in computer readable form, and is hereby 
incorporated by reference herein The computer program 
listing is contained in 6 files, the name, sizes and creation 
date of which are as follows: AUXILARY_FILES.TXT (117K, 
14-NOV-03); BINDING_SITE_SCORING.TXT (17K, 
14-NOV-03); EDIT_DISTANCE.TXT (144K,24-Nov-03); 
FIRST-K.TXT (96K, 24-Nov-03); HAIRPIN_PREDICTION.TXT 
(47K, 14-NOV-03); TWO_PHASED_SIDE_SELECTOR.TXT (4K, 
14-NOV-03); and TWO_PHASED_PREDICTOR.TXT (74K, 

14-NOV-03). 
Summary of Invention 

[0035] The present invention discloses over a thousand novel 

regulatory microRNA (miRNA) oligonucleotides referred to 
here as Genomic Address Messenger (GAM) oligonu- 
cleotides, which GAM oligonucleotides are detectable us- 



ing a novel bioinformatic approach, and go undetected by 
conventional molecular biology methods. Each GAM 
oligonucleotide specifically inhibits translation of one of 
more target genes by hybridization of an RNA transcript 
encoded by the GAM, to a site located in an untranslated 
region (UTR) of the mRNA of one or more of the target 
genes. 

[0036] The present invention represents a scientific break- 
through, disclosing novel mlRNA oligonucleotides the 
number of which is dramatically larger than previously be- 
lieved existed. Prior-art studies reporting mlRNAs ((Lau et 
al.. Science 294:858-862 (2001), Lagos-Quintana et al.. 
Science 294: 853-858 (2001)) discovered 93 miRNAs in 
several species, including 21 in human, using conven- 
tional molecular biology methods, such as cloning and se- 
quencing. Bioinformatic detection of mlRNA oligonu- 
cleotides has not been done prior to the present inven- 
tion: Despite advanced genome projects, computer-as- 
sisted detection of genes encoding functional RNAs re- 
mains problematic (Lagos-Quintana et al.,2001). 

[0037] Molecular biology methodologies employed by these 
studies are limited in their ability to detect rare mlRNA 
oligonucleotides, since these studies relied on sequencing 



of a limited number of clones (300 clones by Lau and 100 
clones by Lagos-Quintana) of small segments (i.e. size 
fractionated) of RNA. MicroRNAs detected in these studies 
therefore, represent the more prevalent among the mlRNA 
oligonucleotide family, and are typically not be much rarer 
than 1% of all small ~20nt-long RNA oligonucleotides 
present in the tissue from the RNA was extracted.. 

[0038] Recent studies state the number of mlRNA genes to be 
limited, and describe the limited sensitivity of available 
methods for detection of miRNA: The estimate of 255 hu- 
man genes is an upper bound implying that no more than 
40 miRNA genes remain to be identified in mammals (Lim 
et al.. Science, 299:1540 (2003)); Estimates place the total 
number of vertebrate miRNA genes at about 200-250 
(Ambros et al. Curr. Biol. 13:807-818 (2003)); and Confir- 
mation of very low abundance mlRNAs awaits the applica- 
tion of detection methods more sensitive than Northern 
blots (Ambros, 2003). 

[0039] The oligonucleotides of the present invention represent a 
revolutionary new dimension of genomics and of biology: 
a dimension comprising a huge number of non-protein 
coding oligonucleotides which modulate expression of 
thousands of proteins and are associated with numerous 



major diseases. This new dimension disclosed by tlie 
present invention dismantles a central dogma that has 
dominated life-sciences during the past 50 years, a 
dogma which has emphasized the importance of protein 
coding regions of the genome, holding non-protein cod- 
ing regions to be of little consequence, often dubbing 
them junk DNA. 

[0040] Indeed, only in recent months has this long held belief as 
to the low importance of non-protein coding regions been 
vocally challenged. As an example, an article titled The 
Unseen Genome - Gems in the Junk (Gibbs,W.W. Sci. Am. 
289:46-53 (2003)) asserts that the failure to recognize 
the importance of non-protein- coding regions may well 
go down as one of the biggest mistakes in the history of 
molecular biology. Gibbs further asserts that what was 
dammned as junk because it was not understood, may in 
fact turn out to be the very basis of human complexity. 
The present invention provides a dramatic leap in under- 
standing specific important roles of non-protein coding 
regions. 

[0041] An additional scientific breakthrough of the present in- 
vention is a novel conceptual model disclosed by the 
present invention, which conceptual model is preferably 



used to encode in a genome the determination of cell- 
differentiation, utilizing oligonucleotides and polynu- 
cleotides of the present invention. 
[0042] In various preferred embodiments, the present invention 
seeks to provide an improved method and system for 
specific modulation of the expression of specific target 
genes involved in significant human diseases.lt also pro- 
vides an improved method and system for detection of the 
expression of novel oligonucleotides of the present inven- 
tion, which modulate these target genes. In many cases 
the target genes may be known and fully characterized, 
however in alternative embodiments of the present inven- 
tion, unknown or less well characterized genes may be 
targeted. 

[0043] Accordingly, the invention provides several substantially 
pure nucleic acids (e.g., genomic DNA, cDNA or synthetic 
DNA) each comprising a novel GAM oligonucleotide, vec- 
tors comprising the DMAs, probes comprising the DMAs, a 
method and system for selectively modulating translation 
of known target genes utilizing the vectors, and a method 
and system utilizing the GAM probes to modulate expres- 
sion of target genes. 

[0044] A Nucleic acid is defined as a ribonucleic acid (RNA) 



molecule, or a deoxyribonucleic acid (DNA) molecule, or 
complementary deoxyribonucleic acid (cDNA), comprising 
either naturally occurring nucleotides or non-naturally oc- 
curring nucleotides. 
[0045] Substantially pure nucleic acid. Isolated Nucleic Acid, Iso- 
lated Oligoucleotide and Isolated Polynucleotide are de- 
fined as a nucleic acid that is free of the genome of the 
organism from which the nucleic acid is derived, and in- 
clude, for example, a recombinant nucleic acid which is 
incorporated into a vector, into an autonomously replicat- 
ing plasmid or virus, or into the genomic nucleic acid of a 
prokaryote or eukaryote at a site other than its natural 
site; or which exists as a separate molecule (e.g., a cDNA 
or a genomic or cDNA fragment produced by PGR or re- 
striction endonuclease digestion) independent of other 
nucleic acids. 

[0046] An Oligonucleotide is defined as a nucleic acid comprising 
2-139 nucleotides, or preferably 16-120 nucleotides. A 
Polynucleotide is defined as a nucleic acid comprising 
140-5000 nucleotides, or preferably 140-1000 nu- 
cleotides. 

[0047] A Complementary sequence is defined as a first nu- 
cleotide sequence which reverse complementary of a sec- 



ond nucleotide sequence: the first nucleotide sequence is 
reversed relative to a second nucleotide sequence, and 
wherein each nucleotide in the first nucleotide sequence is 
complementary to a corresponding nucleotide in the sec- 
ond nucleotide sequence (e.g. ATCGC is the complemen- 
tary sequence of GCCAT). 
[0048] Hybridization, Binding and Annealing are defined as hy- 
bridization, under in-vivo physiologic conditions, of a first 
nucleic acid to a second nucleic acid, which second nu- 
cleic acid is at least partially complementary to the first 
nucleic acid. 

[0049] A Hairpin Sturucture is defined as an oligonucleotide hav- 
ing a nucleotide sequence that is 50-140 nucleotides in 
length, the first half of which nucleotide sequence is at 
least partially complementary to the second part thereof, 
thereby causing the nucleic acid to fold onto itself, form- 
ing a secondary hairpin structure. 

[0050] A Hairpin Shaped Precursor is defined as a Hairpin Struc- 
ture which is processed by a Dicer enzyme complex, 
yielding an oligonucleotide which is about 19 to about 24 
nucleotides in length. 

[0051] "Inhibiting translation" is defined as the ability to prevent 
synthesis of a specific protein encoded by a respective 



gene by means of inhibiting the translation of the mRNA 
of this gene. For example, inhibiting translation may in- 
clude the following steps: (1) a DNA segment encodes an 
RNA, the first half of whose sequence is partially comple- 
mentary to the second half thereof; (2) the precursor folds 
onto itself forming a hairpin-shaped precursor; (3) a Dicer 
enzyme complex cuts the hairpin shaped precursor yield- 
ing an oligonucleotide that isapproximately 22nt in 
length; (4) the oligonucleotide binds complementarily to 
at least one binding site, having a nucleotide sequence 
that is at least partially complentary to the oligonu- 
cleotide, which binding site is located in the mRNA of a 
target gene, preferably in the untranslated region (UTR) of 
a target gene, such that the binding inhibts translation of 
the target protein. 

[0052] A "Translation inhibitor site" is defined as the minimal 
DNA sequence sufficient to inhibit translation. 

[0053] There is thus provided in accordance with a preferred em- 
bodiment of the present invention a bioinformatically de- 
tectable isolated oligonucleotide which is endogenously 
processed from a hairpin-shaped precursor, and anneals 
to a portion of a mRNA transcript of a target gene, 
wherein binding of the oligonucleotide to the mRNA tran- 



script represses expression of tlie target gene, and 
wlierein tlie oligonucleotide has at least 80% sequence 
identity with a nucleotide sequence selected from the 
group consisting of SEQ ID NOs: 1-2686. 

[0054] There is still further provided in accordance with another 
preferred embodiment of the present invention a bioinfor- 
matically detectable isolated oligonucleotide which is en- 
dogenously processed from a hairpin-shaped precursor, 
and anneals to a portion of a mRNA transcript of a target 
gene selected from the group consisting of genes shown 
in Table 11, Row 1, wherein binding of the oligonucleotide 
to the mRNA transcript represses expression of the target 
gene, and wherein the oligonucleotide has at least 80% 
sequence identity with a nucleotide sequence selected 
from the group consisting of SEQ ID NOs: 1-2686. 

[0055] There is additionally provided in accordance with another 
preferred embodiment of the present invention a bioinfor- 
matically detectable isolated oligonucleotide having a nu- 
cleotide sequence selected from the group consisting of 
SEQ ID NOs: 1-2686. 

[0056] There is moreover provided in accordance with another 

preferred embodiment of the present invention a bioinfor- 
matically detectable first oligonucleotide which is a por- 



tion of a mRNA transcript of a target gene, and anneals to 
a second oligonucleotide that is endogenously processed 
from a hairpin precursor, wherein binding of the first 
oligonucleotide to the second oligonucleotide represses 
expression of the target gene, and wherein nucleotide se- 
quence of the second nucleotide is selected from the 
group consisting of SEQ ID NOs: 1-2686. 
[0057] There is further provided in accordance with another pre- 
ferred embodiment of the present invention a bioinfor- 
matically detectable first oligonucleotide which is a por- 
tion of a mRNA transcript of a target gene selected from 
the group consisting of genes shown in Table 11 row 1, 
and anneals to a second oligonucleotide that is endoge- 
nously processed from a hairpin precursor, wherein bind- 
ing of the first oligonucleotide to the second oligonu- 
cleotide represses expression of the target gene, and 
wherein nucleotide sequence of the second nucleotide is 
selected from the group consisting of SEQ ID NOs: 
1-2686. 

[0058] There is still further provided in accordance with another 
preferred embodiment of the present invention a bioinfor- 
matically detectable oligonucleotide having a nucleotide 
sequence selected from the group consisting of SEQ ID 



NOs: 1-2686. 

[0059] There is additionally provided in accordance with another 
preferred embodiment of the present invention a bioinfor- 
matically detectable isolated oligonucleotide which is en- 
dogenously processed from a hairpin-shaped precursor, 
and anneals to a portion of a mRNA transcript of a target 
gene associated with Alzheimers disease, wherein binding 
of the oligonucleotide to the mRNA transcript represses 
expression of the target gene. 

[0060] There is moreover provided in accordance with another 

preferred embodiment of the present invention a bioinfor- 
matically detectable isolated oligonucleotide which is en- 
dogenously processed from a hairpin-shaped precursor, 
and anneals to a portion of a mRNA transcript of a target 
gene that is differentially expressed in a tissue affected by 
Alzheimers disease relative an unaffected tissue, wherein 
binding of the oligonucleotide to the mRNA transcript re- 
presses expression of the target gene. 

[0061] There is further provided in accordance with another pre- 
ferred embodiment of the present invention a bioinfor- 
matically detectable isolated oligonucleotide which is en- 
dogenously processed from a hairpin-shaped precursor, 
and anneals to a portion of a mRNA transcript of a target 



gene, the expression of which target gene correlates with 
Alzheimers disease or susceptibility thereto, wherein 
binding of the oligonucleotide to the mRNA transcript re- 
presses expression of the target gene. 

[0062] There is still further provided in accordance with another 
preferred embodiment of the present invention a bioinfor- 
matically detectable isolated oligonucleotide which is en- 
dogenously processed from a hairpin-shaped precursor, 
and anneals to a portion of a mRNA transcript of a target 
gene responsible for the formation of Alzheimers disease, 
wherein binding of the oligonucleotide to the mRNA tran- 
script represses expression of the target gene. 

[0063] There is additionally provided in accordance with another 
preferred embodiment of the present invention a bioinfor- 
matically detectable isolated oligonucleotide which an- 
neals to a portion of a mRNA transcript of a target gene 
selected from the group consisting of: ACE, ACT, 
ADAMIO, ADAM17, ADAMTS5, APBAl, APBA2, APBB2, 
APLPl, APLP2, APOCl, APP, APPBP2, ARHA, AVP, BACE, 
BACE2, BCHE, BLMH, CASP3, CAVl, CAV2, CHAT, CNTF, 
CRH, CRYAB, CTNND2, CTSB, CTSK, DHCR24, DLST, 
DSCRl, FHL2, FLNB, FOS, GFAP, GLUL, GSK3B, GYPA, 
HADH2, HD, HTR2A, HTR2C, lAPP, IGFl, ILIA, LRPl, 



MAPKIO, MART, MME, MPO, MT3, MTNRIA, NCKAPl, NC- 
STN, NCB, OGDH, OLRl, PLAU, PRND, PSENl, PSEN2, 
RACE, SIOOB, SERPINA3, SLC17A7, SLC6A4, SNCA, 
SNCAIP, TGFBl, TNF, TNFRSF5, TNFRSF6, TRPM2 and UBB, 
wherein binding of tlie oligonucleotide to the mRNA tran- 
script represses expression of the target gene. 
[0064] There is moreover provided in accordance with another 

preferred embodiment of the present invention a bioinfor- 
matically detectable isolated oligonucleotide which is en- 
dogenously processed from a hairpin-shaped precursor, 
and anneals to a portion of a mRNA transcript of a target 
gene selected from the group consisting of: ACE, ACT, 
ADAMIO, ADAM17, ADAMTS5, APBAl, APBA2, APBB2, 
APLPl, APLP2, APOCl, APP, APPBP2, ARHA, AVP, BACE, 
BACE2, BCHE, BLMH, CASP3, CAVl, CAV2, CHAT, CNTF, 
CRH, CRYAB, CTNND2, CTSB, CTSK, DHCR24, DLST, 
DSCRl, FHL2, FLNB, FOS, GFAP, GLUL, GSK3B, GYPA, 
HADH2, HD, HTR2A, HTR2C, lAPP, ICFl, ILIA, LRPl, 
MAPKIO, MAPT, MME, MPO, MT3, MTNRIA, NCKAPl, NC- 
STN, NCB, OGDH, OLRl, PLAU, PRND, PSENl, PSEN2, 
RAGE, SIOOB, SERPINA3, SLC17A7, SLC6A4, SNCA, 
SNCAIP, TCFBl, TNF, TNFRSF5, TNFRSF6, TRPM2 and UBB, 
wherein binding of the oligonucleotide to the mRNA tran- 



script represses expression of the target gene. 
[0065] Tliere is furtlier provided in accordance witli anotlier pre- 
ferred embodiment of the present invention a bioinfor- 
matically detectable isolated oligonucleotide which an- 
neals to a portion of a mRNA transcript of a target gene 
associated with Alzheimers disease, which target gene is 
selected from the group consisting of: ACE, ACT, 
ADAMIO, ADAM17, ADAMTS5, APBAl, APBA2, APBB2, 
APLPl, APLP2, APOCl, APP, APPBP2, ARHA, AVP, BACE, 
BACE2, BCHE, BLMH, CASP3, CAVl, CAV2, CHAT, CNTF, 
CRH, CRYAB, CTNND2, CTSB, CTSK, DHCR24, DLST, 
DSCRl, FHL2, FLNB, FOS, GFAP, GLUL, GSK3B, GYPA, 
HADH2, HD, HTR2A, HTR2C, lAPP, IGFl, ILIA, LRPl, 
MAPKIO, MAPT, MME, MPO, MT3, MTNRIA, NCKAPl, NC- 
STN, NGB, OGDH, OLRl, PLAU, PRND, PSENl, PSEN2, 
RAGE, SIOOB, SERPINA3, SLC17A7, SLC6A4, SNCA, 
SNCAIP, TGFBl, TNF, TNFRSF5, TNFRSF6, TRPM2 and UBB, 
wherein binding of the oligonucleotide to the mRNA tran- 
script represses expression of the target gene, and 
wherein the oligonucleotide has at least 80% sequence 
identity with a nucleotide sequence selected from the 
group consisting of: (a) a sequence selected from the 
group consisting of SEQ ID NOs: 1-36 and 37-2686, and 



(b) the complement of a sequence selected from the 
group consisting of SEQ ID NOs: 2845, 3016, 3025, 3095, 
3164, 3258, 3303, 3304, 3465, 3499, 3547, 3594, 3631, 
4058, 4059, 4261, 4340, 4343, 4384, 4397, 4419, 4448, 
4603, 4616, 4709, 4825, 4843, 4862, 4881, 4896, 5050, 
5325, 5394, 5450, 5458, 5471 and 2687-5643. 
[0066] There is still further provided in accordance with another 
preferred embodiment of the present invention a bioinfor- 
matically detectable isolated oligonucleotide which an- 
neals to a portion of a mRNA transcript of a target gene 
selected from the group consisting of: ACE, ACT, 
ADAMIO, ADAM17, ADAMTS5, APBAl, APBA2, APBB2, 
APLPl, APLP2, APOCl, APP, APPBP2, ARHA, AVP, BACE, 
BACE2, BCHE, BLMH, CASP3, CAVl, CAV2, CHAT, CNTF, 
CRH, CRYAB, CTNND2, CTSB, CTSK, DHCR24, DLST, 
DSCRl, FHL2, FLNB, FOS, GFAP, GLUL, GSK3B, GYPA, 
HADH2, HD, HTR2A, HTR2C, lAPP, IGFl, ILIA, LRPl, 
MAPKIO, MAPT, MME, MPO, MT3, MTNRIA, NCKAPl, NC- 
STN, NGB, OGDH, OLRl, PLAU, PRND, PSENl, PSEN2, 
RAGE, SIOOB, SERPINA3, SLC17A7, SLC6A4, SNCA, 
SNCAIP, TGFBl, TNF, TNFRSF5, TNFRSF6, TRPM2 and UBB, 
wherein binding of the oligonucleotide to the mRNA tran- 
script represses expression of the target gene, and 



wherein the oligonucleotide has at least 80% sequence 
identity with a nucleotide sequence selected from the 
group consisting of SEQ ID NOs: 1-36 and 37-2686. 
[0067] There is additionally provided in accordance with another 
preferred embodiment of the present invention a bioinfor- 
matically detectable isolated oligonucleotide which in- 
cludes from about 19 to about 22 nucleotides, which nu- 
cleotides are partially complementary to a target gene se- 
lected form the group consisting of: ACE, ACT, ADAM 10, 
ADAM17, ADAMTS5, APBAl, APBA2, APBB2, APLPl, APLP2, 
APOCl, APP, APPBP2, ARHA, AVP, BACE, BACE2, BCHE, 
BLMH, CASP3, CAVl, CAV2, CHAT, CNTF, CRH, CRYAB, 
CTNND2, CTSB, CTSK, DHCR24, DLST, DSCRl, FHL2, FLNB, 
FOS, GFAP, GLUL, GSK3B, GYPA, HADH2, HD, HTR2A, 
HTR2C, lAPP, IGFl, ILIA, LRPl, MAPKIO, MAPT, MME, 
MPO, MT3, MTNRIA, NCKAPl, NCSTN, NGB, OGDH, OLRl, 
PLAU, PRND, PSENl, PSEN2, RAGE, SIOOB, SERPINA3, 
SLC17A7, SLC6A4, SNCA, SNCAIP, TGFBl, TNF, TNFRSF5, 
TNFRSF6, TRPM2 and UBB, and wherein the oligonu- 
cleotide is endogenously processed from a hairpin-form 
precursor, and includes at least 19 contiguous nucleotides 
from a sequence selected from the group consisting of 
SEQ ID NOs: 1-36 and 37-2686. 



[0068] There is moreover provided in accordance with another 

preferred embodiment of the present invention a bioinfor- 
matically detectable first oligonucleotide which is a por- 
tion of a mRNA transcript of a target gene selected from 
the group consisting of: ACE, ACT, ADAM 10, ADAM 17, 
ADAMTS5, APBAl, APBA2, APBB2, APLPl, APLP2, APOCl, 
APP, APPBP2, ARHA, AVP, BACE, BACE2, BCHE, BLMH, 
CASP3, CAVl, CAV2, CHAT, CNTF, CRH, CRYAB, CTNND2, 
CTSB, CTSK, DHCR24, DLST, DSCRl, FHL2, FLNB, FOS, 
CFAP, GLUL, GSK3B, GYPA, HADH2, HD, HTR2A, HTR2C, 
lAPP, ICFl, ILIA, LRPl, MAPKIO, MAPT, MME, MPO, MT3, 
MTNRIA, NCKAPl, NCSTN, NGB, OGDH, OLRl, PLAU, 
PRND, PSENl, PSEN2, RAGE, SIOOB, SERPINA3, SLC17A7, 
SLC6A4, SNCA, SNCAIP, TGFBl, TNF, TNFRSF5, TNFRSF6, 
TRPM2 and UBB, and anneals to a second oligonucleotide 
that is endogenously processed from a hairpin precursor, 
wherein binding of the first oligonucleotide to the second 
oligonucleotide represses expression of the target gene. 

[0069] There is further provided in accordance with another pre- 
ferred embodiment of the present invention a bioinfor- 
matically detectable first oligonucleotide which is a por- 
tion of a mRNA transcript of a target gene selected from 
the group consisting of: ACE, ACT, ADAM 10, ADAM 17, 



ADAMTS5, APBAl, APBA2, APBB2, APLPl, APLP2, APOCl, 
APP, APPBP2, ARHA, AVP, BACE, BACE2, BCHE, BLMH, 
CASP3, CAVl, CAV2, CHAT, CNTF, CRH, CRYAB, CTNND2, 
CTSB, CTSK, DHCR24, DLST, DSCRl, FHL2, FLNB, FOS, 
CFAP, CLUL, GSK3B, GYPA, HADH2, HD, HTR2A, HTR2C, 
lAPP, IGFl, ILIA, LRPl, MAPKIO, MAPT, MME, MPO, MT3, 
MTNRIA, NCKAPl, NCSTN, NGB, OGDH, OLRl, PLAU, 
PRND, PSENl, PSEN2, RAGE, SIOOB, SERPINA3, SLC17A7, 
SLC6A4, SNCA, SNCAIP, TGFBl, TNF, TNFRSF5, TNFRSF6, 
TRPM2 and UBB, and anneals to a second oligonucleotide 
that is endogenously processed from a hairpin precursor, 
wherein binding of the first oligonucleotide to the second 
oligonucleotide represses expression of the target gene, 
and wherein nucleotide sequence of the second nucleotide 
is selected from the group consisting of SEQ ID NOs: 1-36 
and 37-2686. 

[0070] There is still further provided in accordance with another 
preferred embodiment of the present invention a bioinfor- 
matically detectable isolated oligonucleotide which an- 
neals to a portion of a mRNA transcript of APP gene, 
wherein binding of the oligonucleotide to the mRNA tran- 
script represses expression of APP gene. 

[0071] There is additionally provided in accordance with another 



preferred embodiment of the present invention a bioinfor- 
matically detectable isolated oligonucleotide which is en- 
dogenously processed from a hairpin-shaped precursor, 
and anneals to a portion of a mRNA transcript of APP 
gene, wherein binding of the oligonucleotide to the mRNA 
transcript represses expression of APP gene. 

[0072] There is moreover provided in accordance with another 

preferred embodiment of the present invention a bioinfor- 
matically detectable isolated oligonucleotide which an- 
neals to a portion of a mRNA transcript of APP gene, 
wherein binding of the oligonucleotide to the mRNA tran- 
script represses expression of APP gene, and wherein the 
oligonucleotide has at least 80% sequence identity with a 
nucleotide sequence selected from the group consisting 
of: (a) a sequence selected from the group consisting of 
SEQ ID NOs: 2,31 and 416-427, and (b) the complement 
of a sequence selected from the group consisting of SEQ 
ID NOs: 3107, 3113 and 3106-3121. 

[0073] There is further provided in accordance with another pre- 
ferred embodiment of the present invention a bioinfor- 
matically detectable isolated oligonucleotide which an- 
neals to a portion of a mRNA transcript of APP gene, 
wherein binding of the oligonucleotide to the mRNA tran- 



script represses expression of APP gene, and wlierein tlie 
oligonucleotide has at least 80% sequence identity with a 
nucleotide sequence selected from the group consisting 
of SEQID NOs: 2, 31 and 416-427. 

[0074] There is still further provided in accordance with another 
preferred embodiment of the present invention a bioinfor- 
matically detectable isolated oligonucleotide which is en- 
dogenously processed from a hairpin-shaped precursor, 
and anneals to a complementary portion of a mRNA tran- 
script of APP gene, wherein binding of the oligonucleotide 
to the mRNA transcript represses expression of APP gene, 
and wherein the oligonucleotide has at least 80% se- 
quence identity with a nucleotide sequence selected from 
the group consisting of SEQ ID NOs: 2, 31 and 416-427. 

[0075] There is additionally provided in accordance with another 
preferred embodiment of the present invention a bioinfor- 
matically detectable isolated oligonucleotide which in- 
cludes from about 19 to about 22 nucleotides, which nu- 
cleotides anneal to a portion of a mRNA transcript of APP 
gene, and wherein the oligonucleotide is endogenously 
processed from a hairpin-form precursor, and includes at 
least 19 contiguous nucleotides from a sequence selected 
from the group consisting of SEQ ID NOs: 2,31 and 



416-427. 

[0076] There is moreover provided in accordance witli another 

preferred embodiment of the present invention a bioinfor- 
matically detectable first oligonucleotide which is a por- 
tion of a mRNA transcript of APP gene, and anneals to a 
second oligonucleotide that is endogenously processed 
from a hairpin precursor, wherein binding of the first 
oligonucleotide to the second oligonucleotide represses 
expression of the target gene. 

[0077] There is further provided in accordance with another pre- 
ferred embodiment of the present invention a bioinfor- 
matically detectable first oligonucleotide which is a por- 
tion of a mRNA transcript of APP gene, and anneals to a 
second oligonucleotide that is endogenously processed 
from a hairpin precursor, wherein binding of the first 
oligonucleotide to the second oligonucleotide represses 
expression of the target gene, and wherein nucleotide se- 
quence of the second nucleotide is selected from the 
group consisting of SEQ ID NOs: 2,31 and 416-427. 

[0078] There is still further provided in accordance with another 
preferred embodiment of the present invention a bioinfor- 
matically detectable isolated polynucleotide which is en- 
dogenously processed into a plurality of hairpin shaped 



precursor oligonucleotides, each of which is endogenously 
processed into a respective oligonucleotide, which in turn 
anneals to a portion of a mRNA transcript of a target 
gene, wherein binding of the oligonucleotide to the mRNA 
transcript represses expression of the target gene. 

[0079] There is additionally provided in accordance with another 
preferred embodiment of the present invention a bioinfor- 
matically detectable isolated oligonucleotide which is en- 
dogenously processed from a hairpin-shaped precursor, 
and anneals to a portion of a mRNA transcript of a target 
gene, wherein binding of the oligonucleotide to the mRNA 
transcript represses expression of the target gene, and 
wherein the target gene does not encode a protein. 

[0080] There is moreover provided in accordance with another 

preferred embodiment of the present invention a bioinfor- 
matically detectable isolated oligonucleotide which is en- 
dogenously processed from a hairpin-shaped precursor, 
and anneals to a portion of a mRNA transcript of a target 
gene, wherein binding of the oligonucleotide to the mRNA 
transcript represses expression of the target gene, and 
wherein a function of the oligonucleotide includes modu- 
lation of cell type. 

[0081] There is further provided in accordance with another pre- 



ferred embodiment of the present invention a bioinfor- 
matically detectable isolated oligonucleotide which is en- 
dogenously processed from a hairpin-shaped precursor, 
and anneals to a portion of a mRNA transcript of a target 
gene, wherein binding of the oligonucleotide to the mRNA 
transcript represses expression of the target gene, and 
wherein the oligonucleotide is maternally transferred by a 
cell to at least one daughter cell of the cell, and a function 
of the oligonucleotide includes modulation of cell type of 
the daughter cell. 
[0082] There is still further provided in accordance with another 
preferred embodiment of the present invention a method 
for bioinformatic detection of microRNA oligonucleotides, 
the method including: bioinformatically detecting a hair- 
pin shaped precursor oligonucleotide, bioinformatically 
detecting an oligonucleotide which is endogenously pro- 
cessed from the hairpin shaped precursor oligonucleotide, 
and bioinformatically detecting a target gene of the 
oligonucleotide wherein the oligonucleotide anneals to at 
least one portion of a mRNA transcript of the target gene, 
and wherein the binding represses expression of the tar- 
get gene, and the target gene is associated with a disease. 
Brief Description of Drawings 



[0083] Fig. 1 is a simplified diagram illustrating a genomic differ- 
entiation enigma that the present invention addresses; 

[0084] Figs. 2, 3 and 4 are schematic diagrams which, when 

taken together, provide an analogy that illustrates a con- 
ceptual model of the present invention, addressing the 
genomic differentiation enigma; 

[0085] Figs. 5A and 5B are schematic diagrams, which when 

taken together, illustrate a "genomic records" concept of 
the conceptual model of the present invention, addressing 
the genomic differentiation enigma; 

[0086] Fig. 6 is a schematic diagram illustrating a 'genomically 
programmed cell differentiation" concept of the concep- 
tual model of the present invention, addressing the ge- 
nomic differentiation enigma; 

[0087] Fig. 7 is a schematic diagram illustrating a "genomically 
programmed cell-specific protein expression modulation" 
concept of the conceptual model of the present invention, 
addressing the genomic differentiation enigma; 

[0088] Fig. 8 is a simplified diagram illustrating a mode by which 
an oligonucleotide of a novel group of oligonucleotides of 
the present invention, modulates expression of known 
target genes; 

[0089] Fig. 9 is a simplified block diagram illustrating a bioinfor- 



matic oligonucleotide detection system capable of detect- 
ing oligonucleotides of the novel group of oligonu- 
cleotides of the present invention, which system is con- 
structed and operative in accordance with a preferred em- 
bodiment of the present invention; 
[0090] Fig. 10 is a simplified flowchart illustrating operation of a 
mechanism for training of a computer system to recog- 
nize the novel oligonucleotides of the present invention, 
which mechanism is constructed and operative in accor- 
dance with a preferred embodiment of the present inven- 
tion; 

[0091] Fig. IIA is a simplified block diagram of a non-coding 
genomic sequence detector constructed and operative in 
accordance with a preferred embodiment of the present 
invention; 

[0092] Fig. IIB is a simplified flowchart illustrating operation of a 
non-coding genomic sequence detector constructed and 
operative in accordance with a preferred embodiment of 
the present invention; 

[0093] Fig. 12A is a simplified block diagram of a hairpin detec- 
tor constructed and operative in accordance with a pre- 
ferred embodiment of the present invention; 

[0094] Fig. 12B is a simplified flowchart illustrating operation of a 



hairpin detector constructed and operative in accordance 
with a preferred embodiment of the present invention; 

[0095] Fig. 13A is a simplified blocl< diagram of a dicer-cut loca- 
tion detector constructed and operative in accordance 
with a preferred embodiment of the present invention; 

[0096] Fig. 13B is a simplified flowchart illustrating training of a 
dicer-cut location detector constructed and operative in 
accordance with a preferred embodiment of the present 
invention; 

[0097] Fig. 13C is a simplified flowchart illustrating operation of 
a dicer-cut location detector constructed and operative in 
accordance with a preferred embodiment of the present 
invention; 

[0098] Fig. 14A is a simplified block diagram of a target gene 

binding-site detector constructed and operative in accor- 
dance with a preferred embodiment of the present inven- 
tion; 

[0099] Fig. 14B is a simplified flowchart illustrating operation of a 
target gene binding-site detector constructed and opera- 
tive in accordance with a preferred embodiment of the 
present invention; 

[0100] Fig. 15 is a simplified flowchart illustrating operation of a 
function & utility analyzer constructed and operative in 



accordance with a preferred embodiment of the present 
invention; 

[0101] Fig. 16 is a simplified diagram describing a novel bioin- 
formatically detected group of regulatory polynucleotides 
referred to here as Genomic Record (GR) polynucleotide, 
each of which encodes an ^operon-like^ cluster of novel 
miRNA-like oligonucleotides, which in turn modulates ex- 
pression of one or more target genes; 

[0102] Fig. 17 is a simplified diagram illustrating a mode by 
which oligonucleotides of a novel group of operon-like 
polynucleotide of the present invention, modulate expres- 
sion of other such polyonucleotides, in a cascading man- 
ner; 

[0103] Fig. 18 is a block diagram illustrating an overview of a 

methodology for finding novel oligonucleotides and novel 
operon-like polynucleotides of the present invention, and 
their respective functions; 

[0104] Fig. 19 is a block diagram illustrating different utilities of 
novel oligonucleotides and novel operon-like polynu- 
cleotides, both of the present invention; 

[0105] Figs. 20A and 20B are simplified diagrams, which when 
taken together illustrate a mode of oligonucleotide-ther- 
apy applicable to novel oligonucleotides of the present in- 



vention; 

[0106] Fig.ZlA is a histogram representing tlie distribution of 
known miRNA oligonucleotides, and that of miRNA-like 
hairpin-shaped oligonucleotides, predicted by the bioin- 
formatics detection engine of the present invention, ex- 
tracted from expressed genome sequences with respect to 
their hairpin detector score. 

[0107] Fig. 2 IB is a table summarizing laboratory validation re- 
sults which validate efficacy of a bioinformatic oligonu- 
cleotide detection system constructed and operative in 
accordance with a preferred embodiment of the present 
invention; 

[0108] Fig. 22A and Fig.22B are a picture and a summary table of 
laboratory results validating the expression of 43 novel 
oligonucleotides detected by a bioinformatic oligonu- 
cleotide detection engine constructed and operative in ac- 
cordance with a preferred embodiment of the present in- 
vention, thereby validating the efficacy of the oligonu- 
cleotide detection engine of the present invention; 

[0109] Fig. 23A is a schematic representation of an "operon-like" 
cluster of novel hairpin sequences detected bioinformati- 
cally by a bioinformatic oligonucleotide detection engine 
constructed and operative in accordance with a preferred 



embodiment of the present invention, and non-GAl\/l hair- 
pin useful as negative controls thereto; 

[01 1 0] Fig. 23B is a schematic representation of secondary fold- 
ing of hairpins of the operon-like cluster of Fig. 23A; 

[0111] Fig. 23C is a picture of laboratory results demonstrating 
expression of novel oligonucleotides of Figs. 23A and 
23B, and lack of expression of the negative controls, 
thereby validating efficacy of bioinformatic detection of 
GAM oligonucleotides and GR polynucleotides of the 
present invention, by a bioinformatic oligonucleotide de- 
tection engine constructed and operative in accordance 
with a preferred embodiment of the present invention; 

[0112] Fig. 24A, is an annotated sequence of EST72223 compris- 
ing known miRNA oligonucleotide MIR98 and novel 
oligonucleotide GAM25 PRECURSOR detected by the 
oligonucleotide detection system of the present invention; 
and 

[0113] Figs. 24B, 24C and 24D are pictures of laboratory results, 
which when taken together demonstrate laboratory con- 
firmation of expression of known oligonucleotide MIR98 
and of novel bioinformatically detected GAM25 RNA re- 
spectively, both of Fig. 24A, thus validating the bioinfor- 
matic oligonucleotide detection system of the present in- 



vention. 
Brief Description of Sequences 



[0114] A Sequence Listing of genomic sequences of tlie present 
invention designated SEQ ID NO: 1 tlirougli SEQ ID NO: 
7,351 is attaclied to tiiis application, and is hereby incor- 
porated herein. The genomic listing comprises the follow- 
ing nucleotide sequences: nucleotide sequences of 1708 
GAMs precursors of respective novel oligonucleotides of 
the present invention; nucleotide sequences of 2686 GAM 
RNA oligonucleotides of respective novel DNA oligonu- 
cleotides of the present invention; and nucleotide se- 
quences of 2957 target gene binding sites of respective 

novel oligonucleotides of the present invention. 
Detailed Description 

[0115] Reference is now made to Fig. 1 which is a simplified dia- 
gram providing a conceptual explanation of a genomic 
differentiation enigma, which the present invention ad- 
dresses, inter alia. 

[0116] Fig. 1 depicts various types of cells in an organism, such 
as a cartilage cell designated by reference numeral 1, a 
liver cell designated by reference numeral 2, a fibroblast 
cell designated by reference numeral 3, and a bone cell 



designated by reference numeral 4, all containing identi- 
cal DNA designated by reference numeral 5. Notwith- 
standing that the various types of cells are all derived 
from an initial fertilized egg cell designated by reference 
numeral 6, each of these cells expresses different proteins 
and accordingly acquires a different shape and function. 

[0117] jhe present invention proposes inter alia that the in- 
evitable conclusion from the foregoing is, however, strik- 
ingly simple: The genome must contain a modular differ- 
entiation coding system. The genome of each cell must 
include multiple modules or records, possibly a different 
one for each cell type, as well as a mechanism causing 
each cell at its inception to be instructed which one of the 
multiple records governs its behavior. 

[0118] This modular code concept may be somewhat difficult to 
grasp, since most persons are accustomed to view things 
from an external viewpoint. An architect, for example, 
looks at a plan of a building, which details exactly where 
each element (block, window, door, electrical switch, etc.) 
is to be placed relative to all other elements, and, using 
the plan, instructs builders to place these elements in 
their designated places. This is an example of an external 
viewpoint: The architect is external to the plan, which it- 



self is external with respect to the physical building, and 
with respect to its various elements. The architect may 
therefore act as an "external organizing agent": seeing the 
full picture and the relationships between all elements, 
and being able to instruct from the outside where to place 
each of them. 

[0119] According to a preferred embodiment of the present in- 
vention, genomic differentiation coding works differently, 
without any such external organizing agent. It comprises 
a smart block (the first cell), which is the architect and the 
plan, and which continuously duplicates itself, somehow 
knowing when to manifest itself as a block and when as a 
window, door, or electrical switch. 

[0120] Reference is now made to Figs. 2A - 4 which are 

schematic diagrams which, when taken together, provide 
an analogy that illustrates a conceptual model of the 
present invention, addressing the genomic differentiation 
enigma. 

[0121] Reference is now made to Fig. 2A. An imaginary talented 
chef, designated by reference numeral 7, is capable of 
preparing any meal provided that he is given specific writ- 
ten cooking instructions. This chef 7 is equipped with two 
items: (a) a thick recipe book, designated by reference nu- 



meral 8, and (b) a small note, designated by reference nu- 
meral 9, having a number scribbled on it. The recipe book 
8 comprises multiple pages, each page detailing how to 
prepare a specific meal. The small note 9 indicates the 
page to be opened, and therefore the meal to be pre- 
pared. The chef looks at the page number written on the 
note, opens the recipe book to the appropriate page, and 
prepares the meal according to the written instructions on 
this page. In the example shown in Fig. 2A, the chef 7 is 
holding a small note 9 bearing the number 127. He there- 
fore opens the book to page 127, as designated by refer- 
ence numeral 10. Since this page contains the recipe for 
preparing bread, the chef 7 prepares a loaf of bread, des- 
ignated by reference numeral 12. Pages of the book, such 
as page 10 in the example shown in Fig. 2A, contains ad- 
ditional information, designated by reference numeral 11 
which additional data is further elaborated hereinbelow 
with reference to Figs. 3 and 4. 
[0122] Reference is now made to Fig. 2B, which depicts two iden- 
tical chefs, a first chef, designated by reference numeral 

13, and a second chef, designated by reference numeral 

14, both holding an identical recipe book, designated by 
reference numeral 8. Although the first chef 13 and the 



second chef 14 are identical, and liold identical recipe 
books 8, they differ in that they hold different small 
notes: the first chef 13 holds a small note designated by 
reference numeral 9, having the number 127 written on it, 
whereas the second chef 14 holds a small note designated 
by reference numeral 15, having the number 134 written 
on it. Accordingly, the first chef 13 opens the book 8 to 
page 127, as designated by reference numeral 10 and, 
based on the instructions written on page 127 prepares a 
loaf of bread, designated by reference numeral 12. The 
second chef 14 opens the book 8 to page 134, as desig- 
nated by reference numeral 16 and, based on the instruc- 
tions written on page 134, prepares a pie, designated by 
reference numeral 17. Pages in the book, such as pages 
10 and 16 in the examples shown in Fig. 2B, contain addi- 
tional information, designated by reference numeral 11 
which additional information is further elaborated herein- 
below with reference to Figs. 3 and 4. 
[0123] Reference is now made to Fig. 3 which illustrates a mode 
by which an imaginary chef can duplicate himself yielding 
two identical chefs, instructing each of the identical dupli- 
cate chefs to prepare a different meal. As an example. Fig. 
3 shows chef 21 duplicating himself, yielding two dupli- 



cate chefs: a first duplicate clief designated by reference 
numeral 22 and a second duplicate chef designated by 
reference numeral 23. The duplicate chefs are identical to 
each other and to chef 21. 

[0124] Like chefs 7 and 13 (Fig. 2A and 2B), Fig. 3 shows chef 21 
holding a recipe book 8 and receiving a note 9 bearing the 
number 127. The chef 21 therefore opens the book 8 to 
page 127, as designated by reference numeral 10, and 
prepares a loaf of bread 12. However, Fig. 3 also elabo- 
rates some of the additional information 11 (Figs. 2A and 
2B) found in page 10: the bottom of page 10, bears two 
numbers, 134 and 157. 

[0125] Chef 21 is trained to perform the following three actions 
when he is finished preparing a meal: (a) Duplicate himself 
yielding two duplicate chefs, the first duplicate chef 22 
and the second duplicate chef 23; (b) Duplicate his recipe 
book 8, handing an identical copy to each of the duplicate 
chefs 22 and 23; and (c) Write down the numbers found at 
the bottom of the page he was instructed to open the 
book to. In the example of chef 21, since he was in- 
structed to open the book to page 10, he writes the num- 
bers 134 and 157 on two respective notes designated by 
reference numerals 15 and 24, and hands note 15 bearing 



the number 134 to the first duplicate chef 22 and note 24 
bearing the number 157 to the second duplicate chef 23. 
[0126] Accordingly, the first duplicate chef 22 receives note 15 
bearing the number 134 and therefore opens the recipe 
book 8 to page 134, as designated by reference numeral 

16, and prepares a pie, designated by reference numeral 

17. The second duplicate chef 23 receives note 24 bearing 
the number 157 and therefore opens the recipe book 8 to 
page 157, as designated by reference numeral 25, and 
prepares rice, designated by reference numeral 26. 

[0127] It is appreciated that while chef 21 and duplicate chefs 22 
and 23 are identical and hold identical recipe books 8, 
they each prepare a different meal. It is also appreciated 
that the meals prepared by the first duplicate chef 22 and 
the second duplicate chef 23 are determined by chef 21, 
and are mediated by the differently numbered notes 15 
and 24 passed on from chef 21 to duplicate chefs 22 and 
23 respectively. 

[0128] It is further appreciated that the mechanism illustrated by 
Fig. 3 enables an unlimited lineage of chefs to divide into 
duplicate, identical chefs and to determine the meals 
those duplicate chefs would prepare. As an example, since 
the first duplicate chef 22 is directed to page 134, as des- 



ignated by reference numeral 16, when he duplicates 
himself (not shown), he will instruct his two duplicate 
chefs to prepare meals specified on pages the numbers of 
which are written at the bottom of page 134, i.e. pages 
114 and 193 respectively. Similarly, the second duplicate 
chef 23 will instruct its duplicate chefs to prepare meals 
specified on pages 121 and 146 respectively, etc. 

[0129] Reference is now made to Fig. 4, which illustrates a mode 
by which a chef can prepare a meal based on instructions 
written in a shorthand format: The main meal-page to 
which a chef is directed by a small note he is given, 
merely contains a list of numbers which further direct him 
to other pages, each specifying how to prepare an ingre- 
dient of that meal. 

[0130] Jo illustrate this shorthand format Fig. 4 shows a chef, 
designated by reference numeral 27, holding the recipe 
book 8 and the note 9 which bears the number 127. The 
chef 27 accordingly opens the recipe book 8 to page 127, 
as designated by reference numeral 10, and based on in- 
structions on this page prepares bread 12. This is similar 
to chefs 7, 13 and 21 of Figs. 2A, 2B and 3 respectively. 

[0131] However, Fig. 4 also further elaborates some of the addi- 
tional information 11 (Figs. 2A and 2B) found in page 10. 



Fig. 4 shows the cooking "instructions" found on page 10 
for mal<ing bread 12 written in a shorthand format, com- 
prising only three numbers, 118, 175 and 183. The chef 
27 writes these numbers on three respective notes desig- 
nated by reference numerals 28 - 30. The notes 28 -30 
are then used to turn to corresponding pages 31 - 33 of 
the book 8, which pages provide instructions for prepara- 
tion of ingredients required for making bread 12: flour 34, 
milk 35 and salt 36. 
[0132] The analogy provided by Figs. 2A - 4 illustrates the con- 
ceptual model of the present invention addressing the ge- 
nomic differentiation enigma, and may be explained as 
follows: The chefs and duplicate chefs 7, 13, 14, 21 - 23 
and 27 (Figs. 2A - 4) in the given analogy represents cells. 
The thick recipe book 8 represents the DNA 5 (Fig. 1). 
Preparing meals such as bread 12, pie 17 or rice 16 (all of 
Fig. 3) represent the cell manifesting itself as a specific 
cell-type, such as cartilage cell 1, liver cell 2, fibroblast 
cell 3, and bone cell 4 (all of Fig. 1). Ingredients of a meal, 
such as flour 34, milk 35 and salt 36 (all of Fig. 4), repre- 
sent proteins typically expressed by a cell of a certain 
cell-type, such as 1 - 4. Like the different chefs of the 
analogy, having the same thick recipe book 8 yet prepar- 



ing different meals, so do different cells in an organism 
contain the same DNA 5 yet manifest themselves as dif- 
ferent cell types, such as 1 - 4, expressing proteins typi- 
cal of these respective cell types. Application of analogy of 
Figs. 2A - 4 to cell-biology is further described hereinbe- 
low with reference to Figs. 5A - 7. 

[0133] Reference is now made to Figs. 5A and 5B which are 

schematic diagrams, which when taken together illustrate 
a Genomic Records concept of the present invention, ad- 
dressing the genomic differentiation enigma. Figs. 5A and 
5B correspond to Figs. 2A and 2B of the chef analogy de- 
scribed hereinabove. 

[0134] An important aspect of the present invention is the Ge- 
nomic Records concept. According to a preferred embodi- 
ment of the present invention the DNA (the thick recipe 
book 8 in the illustration) comprises a very large number 
of Genomic Records (analogous to pages, such as 10, 16 
and 25, in the recipe book) containing the instructions for 
differentiation of a different cell-type, or developmental 
process. Each Genomic Record comprises by a very short 
genomic sequence which functions as a "Genomic Ad- 
dress" of that Genomic Record (analogous to a page num- 
ber, such as the numbers 127, 134 and 157 appearing in 



Fig. 3, in the recipe boolc). At its inception, in addition to 
the DNA, each cell also receives a short RNA segment (the 
scribbled short note, such as 9, 15, 24 of Fig. 3 in the il- 
lustration). This short RNA segment binds complementar- 
ily to a "Genomic Address" sequence of one of the Ge- 
nomic Records, thereby modulating expression of that 
Genomic Record, and accordingly determining the 
cell's-fate (analogous to opening the recipe book 8 to a 
page corresponding to a number on the scribbled note, 
thereby determining the meal to be prepared). A Genomic 
Record may also comprise multiple short RNA segments 
each of which binds complementarily to a target protein 
coding gene, thus modulating expression of this target 
gene (analogous to the shorthand format whereby a page, 
such as 10, points to other pages, such as 31 -33, encod- 
ing various ingredient, such as 34, 35 and 36, all of Fig. 
4). 

[0135] Reference is now made to Fig. 5A. Fig 5A illustrates a cell 
37, having a genome 38. The genome 38 comprises a 
plurality of Genomic Records, some of which Genomic 
Records correlate to specific cell-types. As an example, 6 
such genomic records are shown, corresponding to 6 cell- 
types: LYMPH genomic record 39, FIBROBLAST genomic 



record 40, MUSCLE genomic record 41, BONE genomic 
record 42, CARTILAGE genomic record 43 and NERVE ge- 
nomic record 44. Eacli genomic record comprises genomic 
instructions on differentiation into a specific cell-type, as 
further elaborated hereinbelow with reference to Fig. 7. At 
cell inception, the cell 37 receives a maternal short RNA 
segment 46, which activates one of the genomic records, 
causing the cell to differentiate according to the instruc- 
tions this genomic record comprises. As an example Fig. 
5A illustrates reception of a maternal short RNA segment 
designated 46 having a nucleotide sequence herein sym- 
bolically represented by A' . 
[0136] The FIBROBLAST genomic record 40 contains a binding 
site having a nucleotide sequence symbolically repre- 
sented by A, which is complementary to the nucleotide 
sequence of A' , and therefore the short RNA segment 46 
binds to the FIBROBLAST genomic record 40. This binding 
activates the FIBROBLAST genomic record, causing the cell 
37 to differentiate into a fibroblast cell-type 3 (Fig. 1). 
Other genomic records, designated by reference numerals 
39 and 41 - 44, comprise binding sites having nucleotide 
sequences that are symbolically represented by F, E, B, C 
and D, which are not complementary of the nucleotide se- 



quence of the short RNA segment 46, and are therefore 
not activated thereby. Genomic Records, such as the Fl- 
BR0B1_AST genomic record 40 contain additional informa- 
tion, designated by reference numeral 45, which is further 
elaborated hereinbelow with reference to Figs. 6 and 7. 
[0137] Reference is now made to Fig. 5B, which is a simplified 

schematic diagram, illustrating cellular differentiation me- 
diated by the "Genomic Records" concept. Fig. 5B depicts 
2 cells in an organism, CELL A designated by reference 
numeral 47 and CELL B designated by reference numeral 
48, each having a genome 38. It is appreciated that since 
CELL A 47 and CELL B 48 are cells in the same organism, 
the genome 38 of cells 47 and 48 is identical. Despite 
having an identical genome 38, CELL A 47 differentiates 
differently from CELL B 48, due to activation of different 
genomic records in these two cells. In CELL A 47 the Fl- 
BRO GENOMIC RECORD 40 is activated, causing CELL A 47 
to differentiate into a FIBROBLAST CELL 3, whereas in CELL 
B 48 the BONE GENOMIC RECORD 42 is activated, causing 
the CELL B 48 to differentiate into a BONE CELL 4 (Fig. 1). 
The cause for activation of different genomic records in 
these two cells is the different maternal short RNA which 
they both received: CELL A 47 received a maternal short 



RNA segment designated 46 bearing a nucleotide se- 
quence represented by A' activating genomic record FIBRO 
40, whereas CELL B 48 received a maternal short RNA 
segment designated 49 bearing a nucleotide sequence 
represented by B' activating genomic record BONE 42. 

[0138] Reference is now made to Fig. 6 which is a schematic dia- 
gram illustrating a "genomically programmed cell differ- 
entiation" concept of the conceptual model of the present 
invention, addressing the genomic differentiation enigma. 

[0139] A cell designated CELL A 50 divides into 2 cells desig- 
nated CELL B 51 and CELL C 52. CELL A 50, CELL B 51 and 
CELL C 52 each comprise a GENOME 38, which GENOME 
38 comprises a plurality of GENOMIC RECORDS, herein 
exemplified by reference numerals 40, 42 and 43. It is ap- 
preciated that since CELL A 50, CELL B 51 and CELL C 52 
are cells in the same organism, the GENOME 38 of these 
cells, and the GENOMIC RECORDS, exemplified by 40, 42 
and 43, the genome of these cells comprises, are identical 
in these cells. 

[0140] As described above with reference to Fig. 5B, at its incep- 
tion, CELL A 50 receives a maternal short RNA segment, 
designated by reference numeral 46, having a nucleotide 
sequence represented by A' and outlined by a broken line. 



which activates the FIBRO genomic record 40, thereby 
causing CELL A 50 to differentiate into a FIBROBLAST CELL 
3. However, Fig. 6 elaborates some of the additional in- 
formation 45 (Fig. 5A) of the genomic records: Genomic 
record may also comprise two short genomic sequences, 
referred to here as Daughter Cell Genomic Addresses. 
Blocks designated B and C are Daughter Cell Genomic Ad- 
dresses of the FIBRO Genomic Record. At cell division, 
each parent cell transcribes two short RNA segments, cor- 
responding to the two Daughter Cell Genomic Addresses 
of the Genomic Record of that parent cell, and transfers 
one to each of its two daughter cells. CELL A 50 tran- 
scribes and transfers to its two daughter cells 51 and 52, 
two short RNA segments, designated by reference numer- 
als 49 and 53, outlined by a broken line and designated B' 
and C , corresponding to daughter cell genomic ad- 
dresses designated B and C comprised in the FIBRO ge- 
nomic record 40. 
[0141] CELL B 51 therefore receives the above mentioned mater- 
nal short RNA segment designated 49 having a nucleotide 
sequence represented by B' , which binds complementarily 
to genomic address designated B of the BONE genomic 
record 42, thereby activating this genomic record, which 



in turn causes CELL B 51 to differentiate into a BONE CELL 
4. Similarly, CELL C 52 receives the above mentioned ma- 
ternal short RNA segment designated 53 having a nu- 
cleotide sequence represented by C , which binds com- 
plementarily to genomic address designated C of a CAR- 
TILAGE genomic record 43, thereby activating this ge- 
nomic record, which in turn causes CELL C 52 to differen- 
tiate into a CARTILAGE CELL 1 (Fig. 1). 

[0142] It is appreciated that the mechanism illustrated by Fig. 6 
enables an unlimited lineage of cells to divide into daugh- 
ter cells containing the same DNA 5 (Fig. 1), and to deter- 
mine the cell-fate of these daughter cells. For example, 
when CELL B 51 and CELL C 52 divide into their respective 
daughter cells (not shown), they will transfer short RNA 
segments designated by reference numerals 54 - 57, to 
their respective daughter cells. The cell fate of each of 
these daughter cells is effected by the identity of the ma- 
ternal short RNA segments 54 - 57 they each receive, 
which in turn determine the genomic record activated. 

[0143] Reference is now made to Fig. 7 which is a schematic dia- 
gram illustrating a "genomically programmed cell-specific 
protein expression modulation" concept of the conceptual 
model of the present invention, addressing the genomic 



differentiation enigma. 
[0144] Cell A 58 receives a maternal short RNA segment desig- 
nated 46 having a nucleotide sequence represented by A' , 
which activates the FIBROBLAST genomic record 40, by 
complementarily binding to a binding site this genomic 
record comprises, the nucleotide sequence of which bind- 
ing site is designated A. This is similar to the process 
shown in Fig. 5A. However, Fig. 7 further elaborates some 
of the additional information 45 (Fig. 1). The FIBROBLAST 
genomic record 40 comprises 3 short nucleotide seg- 
ments, having nucleotide sequences symbolically repre- 
sented by 1, 2 and 4 respectively, which encode 3 respec- 
tive short RNA oligonucleotides, designated by reference 
numerals 59 - 61. Each of these short RNA oligonu- 
cleotides modulates expression of a respective one of the 
target genes GENE 1, GENE 2 and GENE 4, designated by 
reference numerals 62 - 64 respectively, by complemen- 
tarily binding to a binding site sequence associated with 
that target gene. In a preferred embodiment of the 
present invention, the modulation of expression of target 
genes such as 62 - 64 comprises translation inhibition of 
target genes by complementarily binding to binding sites 
located in untranslated regions of the target genes. Mod- 



ulation of expression of these genes results in CELL A 58 
differentiating into a FIBROBLAST cell-type 3 (Fig. 1). 
[0^45] It is appreciated that the concept of genomic records each 
comprising a cluster of short RNA segments, which seg- 
ments modulate expression of target genes thereby mod- 
ulating differentiation, is compatible with the clusters of 
miRNA oligonucleotides of the present invention, and 
their translational inhibition of respective target genes by 
means of complementarily binding to binding sites lo- 
cated in the untranslated regions of mRNA of these target 
genes. 

[0146] Reference is now made to Fig. 8, which is a simplified dia- 
gram describing how a plurality of novel bioinformatically 
detectable oligonucleotides of the present invention, re- 
ferred to here as Genomic Address Messenger (GAM) 
oligonucleotides, modulate expression of respective tar- 
get genes. 

[0147] CAM oligonucleotides are novel, bioinformatically de- 
tectable, regulatory, non protein coding, micro RNA 
(miRNA)-like oligonucleotides. The method by which 
GAMs are detected is described hereinbelow with addi- 
tional reference to Figs. 9-15. 

[0148] GAM PRECURSOR DNA is encoded by the human genome. 



CAM TARGET GENE is a human gene encoded by the hu- 
man genome. 

[0149] GAM PRECURSOR DNA encodes a GAM PRECURSOR RNA. 
Similar to miRNA oligonucleotides, CAM PRECURSOR RNA 
does not encode a protein. GAM PRECURSOR RNA folds 
onto itself, forming GAM FOLDED PRECURSOR RNA, which 
has a two-dimensional ^hairpin structure \ As is well 
known in the art, this "hairpin structure \ is typical of by 
mlRNA precursor oligonucleotides, and is due to the fact 
that the nucleotide sequence of the first half of the miRNA 
precursor oligonucleotide is a fully or partially comple- 
mentary sequence of the nucleotide sequence of the sec- 
ond half thereof. By complementary is meant a sequence 
which is reversed and wherein each nucleotide is replaced 
by a complementary nucleotide, as is well known in the art 
(e.g. ATGGC is the complementary sequence of GCCAT). 

[0150] An enzyme complex comprising an enzyme called Dicer 
together with other necessary proteins, herein designated 
as the DICER COMPLEX, "dices" the GAM FOLDED PRE- 
CURSOR RNA yielding a GAM RNA, in the form of a single 
stranded -22 nt long RNA segment. The DICER COMPLEX 
is known in the art to dice a hairpin structured miRNA 
precursor, thereby yielding diced miRNA in the form of a 



short ~22nt RNA segment. 

[0151] GAM TARGET GENE encodes a corresponding messenger 
RNA, designated GAI^ TARGET RNA. GAIVI TARGET RNA 
comprises tliree regions, as is typical of mRNA of a pro- 
tein coding gene: a 5' untranslated region, a protein cod- 
ing region and a 3^ untranslated region, designated 
5^UTR, PROTEIN CODING and 3 UTR respectively. 

[0152] CAM RNA binds complementarily (i.e. hybridizes) to one 
or more target binding sites located in untranslated re- 
gions of GAM TARGET RNA. This complementary binding 
is due to the fact that the nucleotide sequence of GAM 
RNA is a partial or fully complementary sequence of the 
nucleotide sequence of each of the target binding sites. 
As an illustration. Fig. 8 shows three such target binding 
sites, designated BINDING SITE I, BINDING SITE II and 
BINDING SITE III respectively. It is appreciated that the 
number of target binding sites shown in Fig. 8 is only il- 
lustrative and that any suitable number of target binding 
sites may be present. It is further appreciated that al- 
though Fig. 8 shows target binding sites only in the 
3^ UTR region, these target binding sites may be located 
instead in the 5 " UTR region or in both 3 " UTR and 5 " UTR 
regions. 



[0153] jhe complementary binding of CAM RNA to target binding 
sites on CAM TARGET RNA, sucli as BINDING SITE I, BIND- 
ING SITE II and BINDING SITE III, inhibits translation of 
GAM TARGET RNA into GAM TARGET PROTEIN, which is 
shown surrounded by a broken line. 

[0154] It is appreciated that GAM TARGET GENE in fact represents 
a plurality of GAM target genes. The mRNA of each one of 
this plurality of GAM target genes comprises one or more 
target binding sites, each having a nucleotide sequence 
which is at least partly complementary to GAM RNA, and 
which when bound by GAM RNA causes inhibition of 
translation of the GAM target mRNA into a corresponding 
GAM target protein. 

[0155] The mechanism of the translational inhibition exerted by 
GAM RNA on one or more GAM TARGET GENE, may be 
similar or identical to the known mechanism of transla- 
tional inhibition exerted by known miRNA oligonu- 
cleotides. 

[0156] Nucleotide sequences of each of a plurality of GAM 

oligonucleotides described by Fig. 8 and their respective 
genomic sources and chromosomal locations are set forth 
in Tables 1-3, hereby incorporated herein. 

[0157] Nucleotide sequences of GAM PRECURSOR RNAs, and a 



schematic representation of a predicted secondary folding 
of CAM FOLDED PRECURSOR RNAs, of each of a plurality 
of CAM oligonucleotides described by Fig. 8 are set forth 
in Table 4, hereby incorporated herein. 

[0158] Nucleotide sequences of a 'diced ^ CAM RNAof each of a 
plurality of CAM oligonucleotides described by Fig. 8 are 
set forth in Table 5, hereby incorporated herein. 

[0159] Nucleotide sequences of target binding sites, such as 

BINDINC SITE I, BINDINC SITE II and BINDINC SITE III found 
on CAM TARCET RNA, of each of a plurality of CAM 
oligonucleotides described by Fig. 8, and a schematic rep- 
resentation of the complementarity of each of these target 
binding sites to each of a plurality of CAM RNAs described 
by Fig. 8 are set forth in Tables 6 and 7, hereby incorpo- 
rated herein. 

[0160] It is appreciated that specific functions and accordingly 

utilities of each of a plurality of CAM oligonucleotides de- 
scribed by Fig. 8 correlate with, and may be deduced 
from, the identity of the CAM TARCET CENEs inhibited 
thereby, whose . functions are set forth in Table 8, hereby 
incorporated herein. 

[0161] Studies documenting well known correlations between 
each of a plurality of CAM TARCET CENEs of the CAM 



oligonucleotides of Fig. 8, and known functions and dis- 
eases are listed in Table 9, hereby incorporated herein. 

[0162] The present invention discloses a novel group of oligonu- 
cleotides, belonging to the mlRNA-like oligonucleotides 
group, here termed GAM oligonucleotides, for which a 
specific complementary binding has been determined 
bioinformatically. 

[0163] Reference is now made to Fig. 9 which is a simplified 

block diagram illustrating a bioinformatic oligonucleotide 
detection system and method constructed and operative 
in accordance with a preferred embodiment of the present 
invention. 

[0164] An important feature of the present invention is a BIOIN- 
FORMATIC OLIGONUCLEOTIDE DETECTION ENGINE 100, 
which is capable of bioinformatically detecting oligonu- 
cleotides of the present invention. 

[0165] The functionality of the BIOINFORMATIC OLIGONU- 
CLEOTIDE DETECTION ENGINE 100 includes receiving EX- 
PRESSED RNA DATA 102, SEQUENCED DNA DATA 104, and 
PROTEIN FUNCTION DATA 106; performing a complex 
process of analysis of this data as elaborated hereinbelow, 
and based on this analysis provides information, desig- 
nated by reference numeral 108, identifying and describ- 



ing features of novel oligonucleotides. 

[0166] EXPRESSED RNA DATA 102 comprises published ex- 
pressed sequence tags (EST) data, published mRNA data, 
as well as other published RNA data. SEQUENCED DNA 
DATA 104 comprises alphanumeric data representing ge- 
nomic sequences and preferably including annotations 
such as information indicating the location of known pro- 
tein coding regions relative to the genomic sequences. 

[0167] PROTEIN FUNCTION DATA 106 comprises information 

from scientific publications e.g. physiological functions of 
known proteins and their connection, involvement and 
possible utility in treatment and diagnosis of various dis- 
eases. 

[0168] EXPRESSED RNA DATA 102 and SEQUENCED DNA DATA 
104 may preferably be obtained from data published by 
the National Center for Biotechnology Informatiion (NCBI) 
at the National Institute of Health (NIH) QenuthJ.P. (2000). 
Methods Mol. Biol. 132:301-312(2000) , herein incorpo- 
rated by reference). 

[0169] ^ as well as from various other published data sources. 

PROTEIN FUNCTION DATA 106 may preferably be obtained 
from any one of numerous relevant published data 
sources, such as the Online Mendelian Inherited Disease 



In Man (OMIM(TM), Hamosh et al., Nucleic Acids Res. 30: 
52-55(2002)) database developed by John Hopkins Uni- 
versity, and also published by NCBI (2000). 

[0170] Prior to or during actual detection of BIOINFORMATICALLY 
DETECTED GROUP OF NOVEL OLIGONUCLEOTIDES 108 by 
the BIOINFORMATIC OLIGONUCLEOTIDE DETECTION EN- 
GINE 100, BIOINFORMATIC OLIGONUCLEOTIDE DETECTION 
ENGINE TRAINING & VALIDATION FUNCTIONALITY 110 is 
operative. This functionality uses one or more known 
miRNA oligonucleotides as a training set to train the 
BIOINFORMATIC OLIGONUCLEOTIDE DETECTION ENGINE 
100 to bioinformatically recognize mlRNA-like oligonu- 
cleotides, and their respective potential target binding 
sites. BIOINFORMATIC OLIGONUCLEOTIDE DETECTION EN- 
GINE TRAINING &VALIDATION FUNCTIONALITY 110 is fur- 
ther described hereinbelow with reference to Fig. 10. 

[0171] The BIOINFORMATIC OLIGONUCLEOTIDE DETECTION EN- 
GINE 100 preferably comprises several modules which are 
preferably activated sequentially, and are described as 
follows: 

[0172] A NON-CODING GENOMIC SEQUENCE DETECTOR 112 op- 
erative to bioinformatically detect non-protein coding ge- 
nomic sequences. The NON-CODING GENOMIC SEQUENCE 



DETECTOR 112 is further described herein below with ref- 
erence to Figs. IIA and IIB. 

[0173] A HAIRPIN DETECTOR 114 operative to bioinformatically 
detect genomic "hairpin-shaped" sequences, similar to 
CAM FOLDED PRECURSOR RNA (Fig. 8). The HAIRPIN DE- 
TECTOR 114 is further described herein below with refer- 
ence to Figs. 12A and 12B. 

[0174] A DICER-CUT LOCATION DETECTOR 116 operative to 
bioinformatically detect the location on a GAM FOLDED 
PRECURSOR RNA which is enzymatically cut by DICER 
COMPLEX (Fig. 8), yielding diced GAM RNA. The DICER- 
CUT LOCATION DETECTOR 116 is further described herein 
below with reference to Figs. 13A - 13C. 

[0175] A TARGET GENE BINDING-SITE DETECTOR 118 operative 
to bioinformatically detect target genes having binding 
sites, the nucleotide sequence of which is partially com- 
plementary to that of a given genomic sequence, such as 
a nucleotide sequence cut by DICER COMPLEX. The TAR- 
GET GENE BINDING-SITE DETECTOR 118 is further de- 
scribed hereinbelow with reference to Figs. 14A and 14B. 

[0176] A FUNCTION & UTILITY ANALYZER 120 operative to ana- 
lyze function and utility of target genes, in order to iden- 
tify target genes which have a significant clinical function 



and utility. Tlie FUNCTION & UTILITY ANALYZER 120 is 
furtlier described hereinbelow with reference to Fig. 15. 
[0177] According to a preferred embodiment of the present in- 
vention the engine 100 may employ a cluster of 40 PCs ( 
XEON (R) , 2.8GHz, with 80GB storage each), connected by 
Ethernet to 8 servers (2-CPU, XEON (TM) 1.2-2.2GHz, with 
~200GB storage each), combined with an 8-processor 
server (8-CPU, Xeon 550Mhz w/ 8GB RAM) connected via 
2 HBA fiber-channels to an EMC CLARIION (TM) 
100-disks, 3.6 terabyte storage device. A preferred em- 
bodiment of the present invention may also preferably 
comprise software which utilizes a commercial database 
software program, such as MICROSOFT (TM) SQL Server 
2000. It is appreciated that the above mentioned hard- 
ware configuration is not meant to be limiting, and is 
given as an illustration only. The present invention may be 
implemented in a wide variety of hardware and software 
configurations. 

[0178] The present invention discloses 1708 novel oligonu- 
cleotides of the GAM group of oligonucleotides, which 
have been detected bioinformatically, as set forth in Ta- 
bles 1-4, and 246 novel polynucleotides of the GR group 
of polynucleotides , which have been detected bioinfor- 



matically. Laboratory confirmation of 43 bioinformatically 
predicted oligonucleotides of the GAM group of oligonu- 
cleotides , and several bioinformatically predicted polynu- 
cleotides of the CR group of polynucleotides, is described 
hereinbelow with reference to Figs. 21 - 24D. 

[0179] Reference is now made to Fig. 10 which is a simplified 

flowchart illustrating operation of a preferred embodiment 
of the BIOINFORMATIC OLIGONUCLEOTIDE DETECTION 
ENGINE TRAINING &VALIDATION FUNCTIONALITY 110 de- 
scribed hereinabove with reference to Fig. 9. 

[0180] BIOINFORMATIC OLIGONUCLEOTIDE DETECTION ENGINE 
TRAINING & VALIDATION FUNCTIONALITY 110 begins by 
training the BIOINFORMATIC OLIGONUCLEOTIDE DETEC- 
TION ENGINE 100 (Fig. 9) to recognize one or more known 
mlRNA oligonucleotides, as designated by reference nu- 
meral 122. This training step comprises HAIRPIN DETEC- 
TOR TRAINING & VALIDATION FUNCTIONALITY 124, fur- 
ther described hereinbelow with reference to Fig. 12A, 
DICER-CUT LOCATION DETECTOR TRAINING & VALIDA- 
TION FUNCTIONALITY 126, further described hereinbelow 
with reference to Fig. 13A and 13B, and TARGET GENE 
BINDING-SITE DETECTOR TRAINING & VALIDATION FUNC- 
TIONALITY 128, further described hereinbelow with refer- 



ence to Fig. 14A. 

[0181] Next, the BIOINFORMATIC OLIGONUCLEOTIDE DETECTION 
ENGINE TRAINING & VALIDATION FUNCITONALITY 110 is 
operative bioinformatically detect novel oligonucleotides, 
using BIOINFORMATIC OLIGONUCLEOTIDE DETECTION EN- 
GINE 100 (Fig. 9), as designated by reference numeral 130. 
Wet lab experiments are preferably conducted in order to 
validate expression and preferably function of some sam- 
ples of the novel oligonucleotides detected by the BIOIN- 
FORMATIC OLIGONUCLEOTIDE DETECTION ENGINE 100, as 
designated by reference numeral 132. Figs. 22A - 24D il- 
lustrate examples of wet-lab validation of the above men- 
tioned sample novel oligonucleotides bioinformatically 
detected in accordance with a preferred embodiment of 
the present invention. 

[0182] Reference is now made to Fig. IIA which is a simplified 

block diagram of a preferred implementation of the NON- 
CODING GENOMIC SEQUENCE DETECTOR 112 described 
hereinabove with reference to Fig. 9. The NON-PROTEIN 
CODING GENOMIC SEQUENCE DETECTOR 112 preferably 
receives at least two types of published genomic data: EX- 
PRESSED RNA DATA 102 and SEQUENCED DNA DATA 104. 
The EXPRESSED RNA DATA 102 may include, inter alia, EST 



data, EST clusters data, EST genome alignment data and 
mRNA data. Sources for EXPRESSED RNA DATA 102 in- 
clude NCBI dbEST, NCBI UniGene clusters and mapping 
data, and TIGR (Kirkness F. and Kerlavage,A.R., Methods 
Mol. Biol. 69:261-268 (1997)) 

[0183] gene indices. SEQUENCED DNA DATA 104 may include se- 
quence data (FASTA format files), and feature annotations 
(GenBank file format) mainly from NCBI databases. Based 
on the above mentioned input data, the NON-PROTEIN 
CODING GENOMIC SEQUENCE DETECTOR 112 produces a 
plurality of NON-PROTEIN CODING GENOMIC SEQUENCES 
136. Preferred operation of the NON-PROTEIN CODING 
GENOMIC SEQUENCE DETECTOR 112 is described herein- 
below with reference to Fig. IIB 

[0184] Reference is now made to Fig. IIB which is a simplified 
flowchart illustrating a preferred operation of the NON- 
CODING GENOMIC SEQUENCE DETECTOR 112 of Fig. 9. 
Detection of NON-PROTEIN CODING GENOMIC SEQUENCES 
136, generally preferably progresses along one of the fol- 
lowing two paths: 

[0185] A first path for detecting NON-PROTEIN CODING GENOMIC 
SEQUENCES 136 (Fig. IIA) begins with receipt of a plural- 
ity of known RNA sequences, such as EST data. Each RNA 



sequence is first compared witli l<nown protein-coding 
DNA sequences, in order to select only those RNA se- 
quences which are non-protein coding, i.e. intergenic or 
intronic sequences. This can preferably be performed by 
using one of many alignment algorithms known in the art, 
such as BLAST (Altschul et al., J. Mol. Biol. 215:403-410 
(1990)). This sequence comparison preferably also pro- 
vides localization of the RNA sequence on the DNA se- 
quences. 

[0186] Alternatively, selection of non-protein coding RNA se- 
quences and their localization on the DNA sequences can 
be performed by using publicly available EST cluster data 
and genomic mapping databases, such as the UNIGENE 
database published by NCBI or the TIGR database. Such 
databases, map expressed RNA sequences to DNA se- 
quences encoding them, find the correct orientation of 
EST sequences, and indicate mapping of ESTs to protein 
coding DNA regions, as is well known in the art. Public 
databases, such as TIGR, may also be used to map an EST 
to a cluster of ESTs, known in the art as Tentative Human 
Consensus and assumed to be expressed as one seg- 
ment.. Publicly available genome annotation databases, 
such as NCBIs GenBank, may also be used to deduce ex- 



pressed intronic sequences. 

[0187] Optionally, an attempt may be made to "expand" the non- 
protein RNA sequences thus found, by searching for tran- 
scription start and end signals, respectively upstream and 
downstream of the location of the RNA on the DNA, as is 
well known in the art. 

[0188] A second path for detecting NON-PROTEIN CODING GE- 
NOMIC SEQUENCES 136 (Fig. IIA) begins with receipt of 
DNA sequences. The DNA sequences are parsed into non 
protein coding sequences, using published DNA annota- 
tion data, by extracting those DNA sequences which are 
between known protein coding sequences. Next, tran- 
scription start and end signals are sought. If such signals 
are found, and depending on their robustness, probable 
expressed non-protein coding genomic sequences are 
obtained. Such approach is especially useful for identify- 
ing novel GAM oligonucleotides which are found in prox- 
imity to other known mlRNA oligonucleotides, or other 
wet-lab validated GAM oligonucleotides . Since, as de- 
scribed hereinbelow with reference to Fig. 16, GAM 
oligonucleotides are frequently found in clusters, se- 
quences located near known mlRNA oligonucleotides are 
more likely to contain novel GAM oligonucleotides. Op- 



tionally, sequence orthology, i.e. sequence conservation in 
an evolutionary related species, may be used to select ge- 
nomic sequences having a relatively high probability of 
containing expressed novel GAM oligonucleotides. 
[0189] Reference is now made to Fig. 12A which is a simplified 

block diagram of a preferred implementation of the HAIR- 
PIN DETECTOR 114 described hereinabove with reference 
to Fig. 9. 

[0190] The goal of the HAIRPIN DETECTOR 114 is to detect hair- 
pin-shaped genomic sequences, similar to those of known 
mlRNA oligonucleotides. A hairpin- shaped genomic se- 
quence is a genomic sequence, having a first half which is 
at least partially complementary to a second half thereof, 
which causes the halves to folds onto themselves, thereby 
forming a hairpin structure, as mentioned hereinabove 
with reference to Fig. 8. 

[0191] The HAIRPIN DETECTOR 114 (Fig. 9) receives a plurality of 
NON-PROTEIN CODING GENOMIC SEQUENCES 136 (Fig. 
IIA). Following operation of HAIRPIN DETECTOR TRAIN- 
ING & VALIDATION FUNCTIONALITY 124 ( Fig. 10), the 
HAIRPIN DETECTOR 114 is operative to detect and output 
hairpin-shaped sequences, which are found in the NON- 
PROTEIN CODING GENOMIC SEQUENCES 136. The hairpin- 



shaped sequences detected by the HAIRPIN DETECTOR 
114 are designated HAIRPINS STRUCTURES ON GENOMIC 
SEQUENCES 138. A preferred mode of operation of the 
HAIRPIN DETECTOR 114 is described hereinbelow with 
reference to Fig. 12 B. 
[0192] HAIRPIN DETECTOR TRAINING & VALIDATION FUNCTION- 
ALITY 124 includes an iterative process of applying the 
HAIRPIN DETECTOR 114 to known hairpin shaped miRNA 
precursor sequences, calibrating the HAIRPIN DETECTOR 
114 such that it identifies a training set of known hairpin- 
shaped miRNA precursor sequences, as well as other sim- 
ilarly hairpin-shaped sequences. In a preferred embodi- 
ment of the present invention, the HAIRPIN DETECTOR 
TRAINING & VALIDATION FUNCTIONALITY 124 trains the 
HAIRPIN DETECTOR 114 and validates each of the steps of 
operation thereof described hereinbelow with reference to 
Fig. 12B. 

[0193] The HAIRPIN DETECTOR TRAINING & VALIDATION FUNC- 
TIONALITY 124 preferably uses two sets of data: the 
aforesaid training set of known hairpin-shaped miRNA 
precursor sequences, such as hairpin-shaped miRNA pre- 
cursor sequences of 440 miRNA oligonucleotides of H. 
sapiens, M. musculus, C. elegans, C. Brigssae and D. 



Melanogaster, annotated in the RFAM database 
(Griffiths-Jones, 2003), and a large bacl<ground set of 
about 350,000 hairpin-shaped sequences found in ex- 
pressed non-protein coding genomic sequences. The 
bacl<ground set is expected to comprise some valid, pre- 
viously undetected hairpin-shaped miRNA-like precursor 
sequences, and many hairpin-shaped sequences which 
are not hairpin-shaped miRNA-like precursors. 
[0194] In order to validate the performance of the HAIRPIN DE- 
TECTOR 114 (Fig. 9), preferably a variation of the k-fold 
cross validation method (Tom M. Mitchell, Machine Learn- 
ing, McGraw Hill (1997)), is employed. This preferred vali- 
dation method is well adapted to deal with the training 
set, which includes large numbers of similar or identical 
mlRNAs. The training set is therefore preferably initially 
divided into groups of mlRNAs such that any two miRNAs 
that belong to different groups have an Edit Distance 
score of at least D=3, i.e. they differ by at least 3 editing 
steps (Dan Gusfield, Algorithms on strings, trees, and se- 
quences: computer science and computational biology, 
Cambridge University Press, 1997). Next, the groups are 
preferably classified into k sets of groups. Standard k-fold 
cross validation is preferably performed on these sets. 



preferably using k=5, such that the training set and the 
test set include at least one sequence from each of the 
groups. It is appreciated that without the prior grouping, 
standard cross validation methods incorrectly indicate 
much higher performance results for the predictors due to 
the redundancy of training examples within the genome 
of a species and across genomes of different species. 

[0195] In a preferred embodiment of the present invention, using 
the abovementioned validation methodology, the efficacy 
of the HAIRPIN DETECTOR 114 (Fig. 9) is confirmed. For 
example, when a similarity threshold is chosen such that 
90% of the known hairpin-shaped mlRNA precursors are 
successfully predicted, only 11% of the approximately 
342,880 background set of hairpin-shaped sequences are 
predicted to be hairpin-shaped mlRNA-like precursors. 

[0196] Reference is now made to Fig. 12B which is a simplified 
flowchart illustrating preferred operation of the HAIRPIN 
DETECTOR 114 of Fig. 9. The HAIRPIN DETECTOR 114 
preferably initially uses a secondary structure folding al- 
gorithm based on free-energy minimization, such as the 
MFOLD algorithm, described in Mathews et al. J. Mol. Biol. 
288:911-940 (1999) and Zuker, M. Nucleic Acids Res. 31: 
3406-3415. (2003), the disclosure of which is hereby in- 



corporated by reference. This algorithm is operative to 
calculate probable secondary structure folding patterns of 
the NON-PROTEIN CODING GENOMIC SEQUENCES 136 
(Fig. IIA) as well as the free-energy of each of these 
probable secondary folding patterns. The secondary 
structure folding algorithm, such as the MFOLD algorithm 
(Mathews, 1997; Zuker 2003), typically provides a listing 
of the base-pairing of the folded shape, i.e. a listing of 
each pair of connected nucleotides in the sequence. 
[0197] Next, the HAIRPIN DETECTOR 114 analyzes the results of 
the secondary structure folding patterns, in order to de- 
termine the presence and location of hairpin folding 
structures. The goal of this second step is to assess the 
base-pairing listing provided by the secondary structure 
folding algorithm, in order to determine whether the 
base-pairing listing describes one or more hairpin type 
bonding pattern. Preferably, sequence segment corre- 
sponding to a hairpin structure is then separately ana- 
lyzed by the secondary structure folding algorithm in or- 
der to determine its exact folding pattern and free-en- 
ergy. 

[0198] The HAIRPIN DETECTOR 114 then assesses the hairpin 

structures found by the previous step, comparing them to 



hairpin structures of l<nown miRNA precursors, using vari- 
ous cliaracteristic hairpin structure features such as 
length of the hairpin structure, length of the loop of mis- 
matched nucleotides at its center, its free-energy and its 
thermodynamic stability, the amount and type of mis- 
matched nucleotides and the existence of sequence re- 
peat-elements. Only hairpins that bear statistically signifi- 
cant resemblance to the training set of hairpin structures 
of known miRNA precursors, according to the abovemen- 
tioned parameters, are accepted. 
[0199] In a preferred embodiment of the present invention, simi- 
larity to the training set of hairpin structures of known 
miRNA precursors is determined using a "similarity score" 
which is calculated using a weighted sum of terms, where 
each term is a function of one of the abovementioned 
hairpin structure features. The parameters of each func- 
tion are learned from the set of hairpin structures of 
known miRNA precursors, as described hereinabove with 
reference to HAIRPIN DETECTOR TRAINING & VALIDATION 
FUNCTIONALITY 124 (Fig. 10). The weight of each term in 
the similarity score is optimized so as to achieve maxi- 
mized separation between the distribution peaks of simi- 
larity scores validated miRNA-precursor hairpin struc- 



tures, and the distribution of similarity scores of hairpin 
structures detected in the bacl<ground set mentioned 
hereinabove with reference to Fig. 12B. 

[0200] In an alternative preferred embodiment of the present in- 
vention, the step described in the preceding paragraph 
may be split into two stages. A first stage implements a 
simplified scoring method, typically based on threshold- 
ing a subset of the hairpin structure features described 
hereinabove, and may employ a minimum threshold for 
hairpin structure length and a maximum threshold for 
free energy. A second stage is preferably more stringent, 
and preferably employs a full calculation of the weighted 
sum of terms described hereinabove. The second stage 
preferably is performed only on the subset of hairpin 
structures that survived the first stage. 

[0201] The HAIRPIN DETECTOR 114 also attempts to select hair- 
pin structures whose thermodynamic stability is similar to 
that of hairpin structures of known mlRNA precursors. 
This may be achieved in various ways. A preferred em- 
bodiment of the present invention utilizes the following 
methodology, preferably comprising three logical steps: 

[0202] First, the HAIRPIN DETECTOR 114 attempts to group hair- 
pin structures into "families" of closely related hairpin 



structures. As is known in the art, a secondary structure 
folding algorithm typically provides multiple alternative 
folding patterns, for a given genomic sequence and indi- 
cates the free energy of each alternative folding pattern.. 
It is a particular feature of the present invention that the 
HAIRPIN DETECTOR 114 preferably assesses the various 
hairpin structures appearing in the various alternative 
folding patterns and groups hairpin structures which ap- 
pear at identical or similar sequence locations in various 
alternative folding patterns into common sequence loca- 
tion based "families" of hairpins. For example, all hairpin 
structures whose center is within 7 nucleotides of each 
other may be grouped into a family". Hairpin structures 
may also be grouped into a family" if their nucleotide se- 
quences are identical or overlap to a predetermined de- 
gree. 

[0203] It is also a particular feature of the present invention that 
the hairpin structure "families" are assessed in order to 
select only those families which represent hairpin struc- 
tures that are as thermodynamically stable as those of 
hairpin structures of known mlRNA precursors. Preferably 
only families which are represented in at least a selected 
majority of the alternative secondary structure folding 



patterns, typically 65%, 80% or 100% are considered to be 
sufficiently stable. 
[0204] It is an additional particular feature of the present inven- 
tion that the most suitable hairpin structure is selected 
from each selected family. For example, a hairpin struc- 
ture which has the greatest similarity to the hairpin struc- 
tures appearing in alternative folding patterns of the fam- 
ily may be preferred. Alternatively or additionally, the 
hairpin structures having relatively low free energy may be 
preferred. 

[0205] Alternatively or additionally considerations of homology to 
hairpin structures of other organisms and the existence of 
clusters of thermodynamically stable hairpin structures 
located adjacent to each other along a sequence may be 
important in selection of hairpin structures. The tightness 
of the clusters in terms of their location and the occur- 
rence of both homology and clusters may be of signifi- 
cance. 

[0206] Reference is now made to Figs. 13A - 13C which together 
describe the structure and operation of the DICER-CUT 
LOCATION DETECTOR 116, described hereinabove with 
Fig. 9. 

[0207] Fig. 13A is a simplified block diagram of a preferred im- 



plementation of the DICER-CUT LOCATION DETECTOR 
116. The goal of the DICER-CUT LOCATION DETECTOR 
116 is to detect the location in which the DICER COMPLEX, 
described hereinabove with reference to Fig. 8, dices GAM 
FOLDED PRECURSOR RNA, yielding GAM RNA. 

[0208] The DICER-CUT LOCATION DETECTOR 116 therefore re- 
ceives a plurality of HAIRPIN STRUCTURES ON GENOMIC 
SEQUENCES 138 (Fig. 12A), and, following operation of 
DICER-CUT LOCATION DETECTOR TRAINING & VALIDA- 
TION FUNCTIONALITY 126 (Fig. 10), is operative to detect 
a plurality of DICER-CUT SEQUENCES FROM HAIRPIN 
STRUCTURES 140. 

[0209] Reference is now made to Fig. 13B which is a simplified 
flowchart illustrating a preferred implementation of 
DICER-CUT LOCATION DETECTOR TRAINING & VALIDA- 
TION FUNCTIONALITY 126. 

[0210] A general goal of the DICER-CUT LOCATION DETECTOR 
TRAINING & VALIDATION FUNCTIONALITY 126 is to ana- 
lyze the dicer-cut locations of known diced mlRNA on re- 
spective hairpin shaped miRNA precursors in order to de- 
termine a common pattern in these locations, which can 
be used to predict dicer cut locations on GAM folded pre- 
cursor RNAs. 



[0211] jhe dicer-cut locations of l<nown miRNA precursors are 
obtained and studied. Locations of tlie 5' and/or 3' ends 
of the known diced miRNAs are preferably represented by 
their respective distances from the 5' end of the corre- 
sponding hairpin shaped miRNA precursor. Additionally or 
alternatively, the 5' and/or 3' ends of the known diced 
miRNAs are preferably represented by the relationship be- 
tween their locations and the locations of one or more nu- 
cleotides along the hairpin shaped miRNA precursor. Ad- 
ditionally or alternatively, the 5' and/or 3' ends of the 
known diced miRNAs are preferably represented by the 
relationship between their locations and the locations of 
one or more bound nucleotide pairs along the hairpin 
shaped mlRNA precursor. Additionally or alternatively, the 
5' and/or 3' ends of the known diced miRNAs are prefer- 
ably represented by the relationship between their loca- 
tions and the locations of one or more mismatched nu- 
cleotide pairs along the hairpin shaped miRNA precursor. 
Additionally or alternatively, the 5' and/or 3' ends of the 
known diced miRNAs are preferably represented by the 
relationship between their locations and the locations of 
one or more unmatched nucleotides along the hairpin 
shaped miRNA precursor. Additionally or alternatively, lo- 



cations of the 5' and/or 3' ends of the known diced miR- 
NAs are preferably represented by their respective dis- 
tances from the loop located at the center of the corre- 
sponding hairpin shaped miRNA precursor. 

[0212] One or more of the foregoing location metrics may be 

employed in the training and validation functionality. Ad- 
ditionally, metrics related to the nucleotide content of the 
diced mlRNA and/or of the hairpin shaped mlRNA precur- 
sor may be employed. 

[0213] In a preferred embodiment of the present invention, 

DICER-CUT LOCATION DETECTOR TRAINING & VALIDA- 
TION FUNCTIONALITY 126 preferably employs standard 
machine learning techniques known in the art of machine 
learning for analysis of existing patterns in a given "train- 
ing set" of examples. These techniques are capable, to a 
certain degree, of detecting similar patterns in other, pre- 
viously unseen examples. Such machine learning tech- 
niques include, but are not limited to neural networks, 
Bayesian networks. Support Vector Machines (SVM), Ge- 
netic Algorithms, Markovian modeling. Maximum Likeli- 
hood modeling, Nearest Neighbor algorithms, Decision 
trees and other techniques, as is well known in the art. 

[0214] In accordance with one embodiment of the present inven- 



tion, machine learning predictors, sucli as a Support Vec- 
tor IVIachine (SVI\/I) predictor, are applied to the aforemen- 
tioned training set and are operative, for example to test 
every possible nucleotide on a hairpin as a candidate for 
being the 5' end or the 3' end of a diced GAM RNA. More 
preferred machine learning predictors include predictors 
based on Nearest Neighbor, Bayesian modeling, and K- 
nearest-neighbor algorithms. A training set of the known 
miRNA precursor sequences is preferably used for training 
multiple separate classifiers or predictors, each of which 
produces a model for the 5' and/or 3' end locations of a 
diced miRNA with respect to its hairpin precursor. The 
models take into account one or more of the various 
miRNA location metrics described above. 
[0215] Performance of the resulting predictors, evaluated on the 
abovementioned validation set of 440 published miRNAs 
using k-fold cross validation (Mitchell, 1997) with k = 3, 
is found to be as follows: in 70% of known miRNAs 5'-end 
location is correctly determined by an SVM predictor 
within up to 2 nucleotides; a Nearest Neighbor (EDIT DIS- 
TANCE) predictor achieves 56% accuracy (247/440); a 
Two-Phased predictor that uses Bayesian modeling (TWO 
PHASED) achieves 80% accuracy (352/440), when only the 



first phase is used. When the second phase (strand choice) 
is implemented by a nave Bayesian model the accuracy is 
55% (244/440), and when the K-nearest-neighbor model- 
ing is used for the second phase, 374/440 decision are 
made and the accuracy is 65% (242/374). A K-near- 
est-neighbor predictor (FIRST-K) achieves 61% accuracy 
(268/440). The accuracies of all predictors are consider- 
ably higher on top scoring subsets of published mlRNA. 
[0216] Finally, in order to validate the efficacy and accuracy of 
the DICER-CUT LOCATION DETECTOR 116, a sample of 
novel oligonucleotides detected thereby is preferably se- 
lected, and validated by wet lab. Laboratory results vali- 
dating the efficacy of the DICER-CUT LOCATION DETEC- 
TOR 116 are described hereinbelow with reference to Figs. 
21 - 24D. 

[0217] Reference is now made to Fig. 13C which is a simplified 
flowchart illustrating operation of DICER-CUT LOCATION 
DETECTOR 116 (Fig. 9), constructed and operative in ac- 
cordance with a preferred embodiment of the present in- 
vention. The DICER CUT LOCATION DETECTOR 116 
preferably comprises a machine learning computer pro- 
gram module, which is trained to recognize dicer-cut lo- 
cations on known hairpin-shaped miRNA precursors, and 



based on this training, is operable to detect dicer-cut lo- 
cations of novel GAM RNAs (Fig. 8) on GAM FOLDED PRE- 
CURSOR RNAs (Fig. 8). In a preferred embodiment of the 
present invention, the dicer-cut location module prefer- 
ably utilizes machine learning algorithms, such as known 
Support Vector Machine (SVM) and more preferably: 
known Bayesian modeling. Nearest Neighbors, and K- 
nearest-neighbor algorithms. 

[0218] When initially assessing a novel GAM FOLDED PRECURSOR 
RNA, all 19-24 nucleotide long segments thereof are ini- 
tially considered as "potential GAM RNAs", since the dicer- 
cut location is initially unknown. 

[0219] For each such potential GAM RNA, the location of its 5' 
end or the locations of its 5' and 3' ends are scored by at 
least one recognition classifier or predictor. 

[0220] In a preferred embodiment of the present invention, the 
DICER-CUT LOCATION DETECTOR 116 (Fig. 9) may use a 
Support Vector Machine predictor trained on and operat- 
ing on features such as the following: 

[0221] Locations of the 5' and/or 3' ends of the known diced 
mlRNAs, which are preferably represented by their re- 
spective distances from the 5' end of the corresponding 
hairpin shaped miRNA precursor. Additionally or alterna- 



lively, the 5' and/or 3' ends of the known diced miRNAs 
are preferably represented by the relationship between 
their locations and the locations of one or more nu- 
cleotides along the hairpin shaped mlRNA precursor. Ad- 
ditionally or alternatively, the 5' and/or 3' ends of the 
known diced miRNAs are preferably represented by the 
relationship between their locations and the locations of 
one or more bound nucleotide pairs along the hairpin 
shaped mlRNA precursor. Additionally or alternatively, the 
5' and/or 3' ends of the known diced miRNAs are prefer- 
ably represented by the relationship between their loca- 
tions and the locations of one or more mismatched nu- 
cleotide pairs along the hairpin shaped miRNA precursor. 
Additionally or alternatively, the 5' and/or 3' ends of the 
known diced miRNAs are preferably represented by the 
relationship between their locations and the locations of 
one or more unmatched nucleotides along the hairpin 
shaped mlRNA precursor. Additionally or alternatively, lo- 
cations of the 5' and/or 3' ends of the known diced miR- 
NAs are preferably represented by their respective dis- 
tances from the loop located at the center of the corre- 
sponding hairpin shaped miRNA precursor; and secondar- 
ily 



[0222] Metrics related to the nucleotide content of the diced 
miRNA and/or of the hairpin shaped miRNA precursor. 

[0223] In another preferred embodiment of the present inven- 
tion, the DICER-CUT LOCATION DETECTOR 116 (Fig. 9) 
preferably employs an "EDIT DISTANCE" predictor, which 
seeks sequences that are similar to those of known miR- 
NAs, utilizing a Nearest Neighbor algorithm, where a sim- 
ilarity metric between two sequences is a variant of the 
Edit Distance algorithm (Dan Cusfield, Algorithms on 
strings, trees, and sequences: computer science and com- 
putational biology, Cambridge University Press, 1997). 
The EDIT DISTANCE predictor is based on an observation 
that miRNA oligonucleotides tend to form clusters , the 
members of which show marked sequence similarity. 

[0224] In yet another preferred embodiment of the present in- 
vention, the DICER-CUT LOCATION DETECTOR 116 (Fig. 9) 
preferably uses a "TWO PHASE" predictor, which predicts 
the dicer-cut location in two distinct phases: (a) selecting 
a double-stranded segment of the GAM FOLDED PRECUR- 
SOR RNA (Fig. 8) comprising the GAM RNA by nave 
Bayesian modeling and (b) detecting which strand of the 
double-stranded segment contains GAM RNA (Fig. 8) by 
employing either nave or by K-nearest-neighbor model- 



ing. K-nearest-neighbor modeling is a variant of tlie 
'FIRST-K' predictor described liereinbelow, witli parame- 
ters optimized for tiiis specific task. The 'TWO PHASE' pre- 
dictor may be operated in two modes: either utilizing only 
the first phase and thereby producing two alternative 
dicer-cut location predictions, or utilizing both phases 
and thereby producing only one final dicer-cut location. 

[0225] In still another preferred embodiment of the present in- 
vention, the DICER-CUT LOCATION DETECTOR 116 
preferably uses a "FIRST-K" predictor, which utilizes a K- 
nearest-neighbor algorithm. The similarity metric between 
any two sequences is 1- E/L , where Lisa parameter, 
preferably 8-10 and E is the edit distance between the two 
sequences, taking into account only the first L nucleotides 
of each sequence. If the K-nearest-neighbor scores of two 
or more locations on the GAM FOLDED PRECURSOR RNA 
(Fig. 8) are not significantly different, these locations are 
further ranked by a Bayesian model, similar to the one de- 
scribed hereinabove. 

[0226] The TWO PHASE and FIRST-K predictors preferably are 

trained on and operate on features such as the following: 

[0227] Locations of the 5' and/or 3' ends of the the known diced 
miRNAs, which are preferably represented by their re- 



spective distances from the 5' end of the corresponding 
hairpin shaped miRNA precursor. Additionally or alterna- 
tively, the 5' and/or 3' ends of the known diced mlRNAs 
are preferably represented by the relationship between 
their locations and the locations of one or more nu- 
cleotides along the hairpin shaped miRNA precursor. Ad- 
ditionally or alternatively, the 5' and/or 3' ends of the 
known diced mlRNAs are preferably represented by the 
relationship between their locations and the locations of 
one or more bound nucleotide pairs along the hairpin 
shaped miRNA precursor. Additionally or alternatively, the 
5' and/or 3' ends of the known diced mlRNAs are prefer- 
ably represented by the relationship between their loca- 
tions and the locations of one or more mismatched nu- 
cleotide pairs along the hairpin shaped miRNA precursor. 
Additionally or alternatively, the 5' and/or 3' ends of the 
known diced miRNAs are preferably represented by the 
relationship between their locations and the locations of 
one or more unmatched nucleotides along the hairpin 
shaped miRNA precursor. Additionally or alternatively, lo- 
cations of the 5' and/or 3' ends of the the known diced 
miRNAs are preferably represented by their respective 
distances from the loop located at the center of the corre- 



spending hairpin sliaped miRNA precursor; and secondar- 
ily 

[0228] Metrics related to the nucleotide content of the diced 
mlRNA and/or of the hairpin shaped miRNA precursor. 

[0229] In accordance with an embodiment of the present inven- 
tion scores of two or more of the abovementioned classi- 
fiers or predictors are integrated, yielding an integrated 
score for each "potential GAM RNA". As an example, Fig. 
13C illustrates integration of scores from two classifiers, a 
3' end recognition classifier and a 5' end recognition clas- 
sifier, the scores of which are integrated to yield an inte- 
grated score. Most preferably, the INTEGRATED SCORE of 
Fig. 13C preferably implements a "best-of-breed" ap- 
proach employing a pair of classifiers and accepting only 
"potential GAM RNAs" that score highly on one of the 
above mentioned "EDIT DISTANCE", or "TWO-PHASE" pre- 
dictors. In this context, "high scores" means scores which 
have been demonstrated to have low false positive value 
when scoring known miRNA oligonucleotides. Alterna- 
tively, the INTEGRATED SCORE may be derived from oper- 
ation of more or less than two classifiers. 

[0230] The INTEGRATED SCORE is evaluated as follows: (a) the 
"potential GAM RNA" having the highest score is prefer- 



ably taken to be the most probable GAM RNA, and (b) if 
the integrated score of this most probable GAM RNA is 
higher than a pre-defined threshold, then the most prob- 
able GAM RNA is accepted as a PREDICTED GAM RNA. 
Preferably, this evaluation technique is not limited to the 
highest scoring potential GAM RNA. 

[0231] Reference is now made to Fig. 14A which is a simplified 
block diagram of a preferred implementation of the TAR- 
GET GENE BINDING-SITE DETECTOR 118 described herein- 
above with reference to Fig. 9. The goal of the TARGET 
GENE BINDING-SITE DETECTOR 118 is to detect one or 
more binding sites such as BINDING SITE I, BINDING SITE II 
and BINDING SITE III (Fig. 8) located in untranslated re- 
gions of the mRNA of a known gene, the nucleotide se- 
quence of which binding sites is partially or fully comple- 
mentary to a GAM RNA, thereby determining that the 
above mentioned known gene is a target gene thereof. 

[0232] The TARGET GENE BINDING-SITE DETECTOR 118 (Fig. 9) 

receives a plurality of DICER-CUT SEQUENCES FROM HAIR- 
PIN STRUCTURES 140 (Fig. 13A), and a plurality of POTEN- 
TIAL TARGET GENE SEQUENCES 142 which are derived 
from SEQUENCED DNA DATA 104 (Fig. 9). 

[0233] TARGET GENE BINDING-SITE DETECTOR TRAINING & VALI- 



DATION FUNCTIONALITY 128 (Fig. 10) is operative to train 
the TARGET GENE BINDING-SITE DETECTOR on known 
mlRNAs and tlieir respective target genes. A sequence 
comparison of sequences of l<nown miRNA oligonu- 
cleotides to sequences of known binding sites of known 
target thereof is performed by utilizing BLAST or other al- 
gorithms such as EDIT DISTANCE. 

[0234] jhe results are preferably employed to define a threshold 
based on scoring distinctions between known miRNA 
binding sites and sequences which are known not to be 
miRNA binding sites. This threshold is used during opera- 
tion of TARGET GENE BINDING-SITE DETECTOR 118 to 
distinguish between mlRNA-like binding sites of potential 
GAM RNA and other sequences. 

[0235] Next, the binding sites are expanded, and determinations 
are made whether if nucleotide sequences immediately 
adjacent to the binding sites found by the sequence com- 
parison algorithm (e.g. BLAST or EDIT DISTANCE), may im- 
prove the match. Free-energy and spatial structure are 
computed for the resulting binding sites. Binding sites 
which are clustered are strongly preferred and binding 
sites found in evolutionarily conserved sequences may 
also be preferred. Free energy, spatial structure and the 



above preferences are reflected in scoring. 
[0236] The resulting scores, characteristic of known binding sites 
(e.g. binding sites of l<nown miRNA oligonucleotides Un-4 
and Let-7 to target genes Lin-14, Lin-41, Lin 28 etc.), 
may be employed for detection of binding-sites of novel 
GAM RNAs. 

[0237] Following operation of TARGET GENE BINDING-SITE DE- 
TECTOR TRAINING & VALIDATION FUNCTIONALITY 128 
(Fig. 10), the TARGET GENE BINDING-SITE DETECTOR 118 
is operative to detect a plurality of POTENTIAL NOVEL 
TARGET GENES HAVING BINDING-SITE/S 144 the nu- 
cleotide sequence of which is partially or fully comple- 
mentary to that of each of the plurality of DICER-CUT SE- 
QUENCES FROM HAIRPIN STRUCTURES 140. Preferred op- 
eration of the TARGET GENE BINDING-SITE DETECTOR 118 
is further described hereinbelow with reference to Fig. 
14B. 

[0238] Reference is now made to Fig. 14B which is a simplified 
flowchart illustrating a preferred operation of the TARGET 
GENE BINDING-SITE DETECTOR 118 of Fig. 9. In a pre- 
ferred embodiment of the present invention, the TARGET 
GENE BINDING-SITE DETECTOR 118 employs a sequence 
comparison algorithm such as BLAST in order to compare 



the nucleotide sequence of each of the plurality of DICER- 
CUT SEQUENCES FROM HAIRPIN STRUCTURES 140 (Fig. 
13A), to the POTENTIAL TARGET GENE SEQUENCES 142 
(Fig.l4A), such as untranslated regions of known mRNAs, 
in order to find crude potential matches. Alternatively, the 
sequence comparison may be performed using a se- 
quence match search tool that is essentially a variant of 
the EDIT DISTANCE algorithm described hereinabove with 
reference to Fig. 13C, and the Nearest Neighbor algo- 
rithm. 

[0239] A sequence comparison of DICER-CUT SEQUENCES FROM 
HAIRPIN STRUCTURES 140 (Fig. 14A) are compared to PO- 
TENTIAL TARGET GENE SEQUENCES 142 (Fig. 14A) by uti- 
lizing BLAST or other algorithms such as EDIT DISTANCE. 

[0240] The results are preferably filtered according to a threshold 
determined in accordance with the scoring resulting from 
the sequence comparison carried out by the TARGET GENE 
BINDING-SITE DETECTOR TRAINING & VALIDATION FUNC- 
TIONALITY 128. 

[0241] Next the binding sites are expanded, and determinations 
are made whether if nucleotide sequences immediately 
adjacent to the binding sites found by the sequence com- 
parison algorithm (e.g. BLAST or EDIT DISTANCE), may im- 



prove the match. 

[0242] Free-energy and spatial structure are computed for the 
resulting binding sites. Binding sites which are clustered 
are strongly preferred and binding sites found in evolu- 
tionarily conserved sequences may also be preferred. Free 
energy, spatial structure and the above preferences are 
reflected in scoring. 

[0243] The resulting scores are compared with scores character- 
stic of known binding sites (e.g. binding sites of known 
miRNA oligonucleotides Lin-4 and Let-7 to target genes 
Lin-14, Lin-41, Lin 28 etc.). 

[0244] For each candidate binding site a score, here termed 

Binding Site Prediction Accuracy, is calculated which esti- 
mates its similarity to known binding sites. This score is 
based on GAM binding site characteristics including, but 
not limited to: 

[0245] The free energy of binding of the GAM RNA - GAM RNA 
binding site complex; 

[0246] Additionally or alternatively, the 5' and/or 3' ends of the 
GAM RNA, preferably represented by the relationship be- 
tween their locations and the locations of one or more nu- 
cleotides along the GAM RNA; Additionally or alternatively, 
the 5' and/or 3' ends of the GAM RNA, preferably repre- 



sented by the relationship between their locations and the 
locations of one or more bound nucleotide pairs along the 
GAM RNA binding site complex; Additionally or alterna- 
tively, the 5' and/or 3' ends of the GAM RNA, preferably 
represented by the relationship between their locations 
and the locations of one or more mismatched nucleotide 
pairs along the GAM RNA binding-site complex; Addition- 
ally or alternatively, the 5' and/or 3' ends of the GAM RNA, 
preferably represented by the relationship between their 
locations and the locations of one or more unmatched nu- 
cleotides along the GAM RNA binding-site complex. 
[0247] In accordance with another preferred embodiment of the 
present invention, binding sites are searched by a re- 
versed process. Sequences of K (preferably 22) nu- 
cleotides of a untranslated regions of a target gene are 
assessed as potential binding sites. A sequence compari- 
son algorithm, such as BLAST or EDIT DISTANCE, is then 
used to search elsewhere in the genome for partially or 
fully complementary sequences which are found in known 
miRNA oligonucleotides or computationally predicted GAM 
oligonucleotides. Only complementary sequences, which 
meet predetermined spatial structure and free energy cri- 
teria as described hereinabove are accepted. Clustered 



binding sites are strongly preferred and potential binding 
sites and potential GAM oligonucleotides which occur in 
evolutionarily conserved genomic sequences are also pre- 
ferred. Scoring of candidate binding sites takes into ac- 
count free energy and spatial structure of the binding site 
complexes, as well as the aforesaid preferences. 

[0248] Target binding sites identified by the TARGET GENE BIND- 
ING-SITE DETECTOR 118 (Fig. 9), are preferably divided 
into four groups: 

[0249] a) binding sites which are exactly complementary to the 
predicted GAM RNA. (Int. mismatch is allowed) 

[0250] b) binding sites which are not exactly complementary to 
the predicted GAM RNA and having 0.8=< Binding Site 
Prediction Accuracy < 1; 

[0251] c) binding sites which are not exactly complementary to 
the predicted GAM RNA and having 0.7= < Binding Site 
Prediction Accuracy < 0.8; and 

[0252] d) binding sites which are not exactly complementary to 
the predicted GAM RNA and having 0.6= < Binding Site 
Prediction Accuracy < 0.7. 

[0253] The average number of mismatched nucleotides in the 
alignment of predicted GAM RNA and a corresponding 
target gene binding-site is smallest in category a and 



largest in category d. 
[0254] In accordance with a preferred embodiment of tlie present 
invention tliere is provided a binding site specific ranking, 
indicative of tlie degree of similarity of characteristics of 
the binding of a GAM to a target gene binding site, to 
binding characteristic of known miRNAs. This ranking 
preferably utilizes the evaluation criteria described here- 
inabove. 

[0255] In accordance with another preferred embodiment of the 
present invention, there is provided a UTR specific ranking 
of CAM to target gene binding. , indicative of the degree 
of similarity of characteristics of the binding of a GAM to a 
cluster of target gene binding sites on a UTR, to binding 
characteristics of known miRNAs to UTRs of correspond- 
ing mlRNA target genes. This ranking preferably is a 
weighted sum of the binding site specific rankings of vari- 
ous clustered binding sites. 

[0256] Reference is now made to Fig. 15 which is a simplified 
flowchart illustrating a preferred operation of the FUNC- 
TION & UTILITY ANALYZER 120 described hereinabove 
with reference to Fig. 9. The goal of the FUNCTION & 
UTILITY ANALYZER 120 is to determine if a potential tar- 
get gene is in fact a valid clinically useful target gene. 



Since a potential novel GAM oligonucleotide binding a 
binding site in the UTR of a target gene is understood to 
inhibit expression of that target gene, and if that target 
gene is shown to have a valid clinical utility, then in such a 
case it follows that the potential novel oligonucleotide it- 
self also has a valid useful function which is the opposite 
of that of the target gene. 
[0257] The FUNCTION & UTILITY ANALYZER 120 preferably re- 
ceives as input a plurality of POTENTIAL NOVEL TARGET 
GENES HAVING BINDING-SITE/S 144 (Fig. 14A), generated 
by the TARGET GENE BINDING-SITE DETECTOR 118 (Fig. 
9). Each potential oligonucleotide is evaluated as follows: 
First, the system checks to see if the function of the po- 
tential target gene is scientifically well established. Prefer- 
ably, this can be achieved bioinformatically by searching 
various published data sources presenting information on 
known function of proteins. Many such data sources exist 
and are published as is well known in the art. Next, for 
those target genes the function of which is scientifically 
known and is well documented, the system then checks if 
scientific research data exists which links them to known 
diseases. For example, a preferred embodiment of the 
present invention utilizes the OMIM(TM) (Hamosh et al. 



2002) database published by NCBI, which summarizes re- 
search publications relating to genes which have been 
shown to be associated with diseases. Finally, the specific 
possible utility of the target gene is evaluated. While this 
process too may be facilitated by bioinformatic means, it 
might require manual evaluation of published scientific 
research regarding the target gene, in order to determine 
the utility of the target gene to the diagnosis and or treat- 
ment of specific disease. Only potential novel oligonu- 
cleotides, the target genes of which have passed all three 
examinations, are accepted as novel oligonucleotide. 
[0258] Reference is now made to Fig. 16, which is a simplified di- 
agram describing each of a plurality of novel bioinformat- 
ically detected regulatory polynucleotide, referred to here 
as Genomic Record (GR) polynucleotide which encodes an 
"operon-like^ cluster of novel micro RNA-like oligonu- 
cleotides each of which in turn modulates expression of at 
least one target gene, the function and utility of which at 
least one target gene is known in the art. GR GPREGURSOR 
DNA is a novel bioinformatically detected regulatory, non 
protein coding, polynucleotide. The method by which GR 
polynucleotide as detected is described hereinabove with 
additional reference to Figs. 9-18. GR GPREGURSOR DNA 



encodes GR PRECURSOR RNA, an RNA molecule, typically 
several hundreds to several thousands nucleotides long. 
GR PRECURSOR RNA folds spatially, forming GR FOLDED 
PRECURSOR RNA. It is appreciated that GR FOLDED PRE- 
CURSOR RNA comprises a plurality of what is known in the 
art as ^ hairpin^ structures. These ^ hairpin^ structures are 
due to the fact that the nucleotide sequence of GR PRE- 
CURSOR RNA comprises a plurality of segments, the first 
half of each such segment having a nucleotide sequence 
which is at least a partial or accurate complementary se- 
quence of the second half thereof, as is well known in the 
art. GR FOLDED PRECURSOR RNA is naturally processed by 
cellular enzymatic activity into separate GAM precursor 
RNAs, herein schematically represented by GAMl FOLDED 
PRECURSOR RNA through GAM3 FOLDED PRECURSOR RNA, 
each of which GAM precursor RNAs being a hairpin 
shaped RNA segment, corresponding to GAM FOLDED 
PRECURSOR RNA of Fig. 8. The above mentioned GAM pre- 
cursor RNAs are diced by DICER COMPLEX of Fig. 8, yield- 
ing respective short RNA segments of about 22 nu- 
cleotides in length, schematically represented by GAMl 
RNA through GAM3 RNA, each of which GAM RNAs corre- 
sponding to GAM RNA of Fig. 8. GAMl RNA, GAM2 RNA 



and GAM3 RNA, each bind complementarily to binding 
sites located in untranslated regions of respective target 
genes, designated CAMl-TARGET RNA, CAM2-TARCET 
RNA and GAM3-TARGET RNA, respectively, which target 
binding site corresponds to a target binding site such as 
BINDING SITE I, BINDING SITE II or BINDING SITE III of 
Fig. 8. This binding inhibits translation of the respective 
target proteins designated GAMl-TARGET PROTEIN, 
GAM2-TARGET PROTEIN and GAM3-TARGET PROTEIN re- 
spectively. It is appreciated that specific functions, and 
accordingly utilities, of each GR polynucleotides of the 
present invention, correlates with, and may be deduced 
from, the identity of the target genes, which are inhibited 
by GAM RNAs comprised in the ^operon-like" cluster of 
said GR polynucleotide schematically represented by 
GAMl TARGET PROTEIN through GAMS TARGET PROTEIN. 
[0259] A listing of GAM oligonucleotide comprised in each of a 

plurality of GR polynucleotide of Fig. 16 is provided in Ta- 
ble 10, hereby incorporated herein. Nucleotide sequences 
of each said GAM oligonucleotide and their respective ge- 
nomic source and chromosomal location are further de- 
scribed hereinbelow with reference to Table 3 hereby in- 
corporated herein. GAM TARGET GENEs of each of said 



GAM oligonucleotides are elaborated hereinbelow with 
reference to Table 7, hereby incorporated herein. The 
functions of each of said CAM TARGET GENEs and their 
association with various diseases, and accordingly the 
utilities of said each of GAM oligonucleotides and hence 
the functions and utilities of each of said GR polynu- 
cleotides are elaborated hereinbelow with reference to Ta- 
ble 8, hereby incorporated herein. Studies establishing 
known functions of each of said GAM TARGET GENEs, and 
correlation of each of said GAM TARGET GENEs to known 
diseases are listed in Table 9, and are hereby incorporated 
herein. 

[0260] The present invention discloses 246 novel genes of the GR 
group of polynucleotides, which have been detected 
bioinformatically. Laboratory confirmation of 2 polynu- 
cleotides of the GR group of polynucleotides is described 
hereinbelow with reference to Figs. 23A - 24D. 

[0261] In summary, the current invention discloses a very large 
number of novel GR polynucleotides each of which en- 
codes a plurality of GAM oligonucleotides, which in turn 
may modulate expression of a plurality of target proteins. 
It is appreciated therefore that the function of GR polynu- 
cleotides is in fact similar to that of the Genomic Records 



concept of the present invention addressing the differen- 
tiation enigma, described hereinabove with reference to 
Fig. 7. 

[0262] Reference is now made to Fig. 17 which is a simplified di- 
agram illustrating a mode by which oligonucleotides of a 
novel group of operon-like polynucleotide described 
hereinabove with reference to Fig. 16 of the present in- 
vention, modulate expression of other such polynu- 
cleotide, in a cascading manner.GRl PRECURSOR DNA and 
GR2 PRECURSOR DNA are two polynucleotides of the novel 
group of operon-like polynucleotides designated GR PRE- 
CURSOR DNA(Fig. 16). As is typical of polynucleotides of 
the GR group of polynucleotides GRl PRECURSOR DNA 
and GR2 PRECURSOR DNA, each encode a long RNA pre- 
cursor, which in turn folds into a folded RNA precursor 
comprising multiple hairpin shapes, and is cut into re- 
spective separate hairpin shaped RNA segments, each of 
which RNA segments being diced to yield a n oligonu- 
cleotide of a group of oligonucleotide designated GAM 
RNA. In this manner GRl yields GAMl RNA, GAM2 RNA 
and GAM3 RNA, and GR2 yields GAM4 RNA, GAMS RNA 
and GAM6 RNA. As Fig. 17 shows, GAMS RNA , which de- 
rives from GRl, binds a binding site located adjacent to 



GR2 GPRECURSOR DNA thus modulating expression of 
GR2, tliereby invol<ing expression of GAI\/I4 RNA , GAIV15 
RNA and GAM6 RNA wliicli derive from GR2. It is appreci- 
ated that the mode of modulation of expression presented 
by Fig. 17 enables an unlimited 'cascading effect' in 
which a GR polynucleotide comprises multiple GAM 
oligonucleotides each of which may modulate expression 
of other GR polynucleotides each such GR polynucleotides 
comprising additional GAM oligonucleotide etc., whereby 
eventually certain GAM oligonucleotides modulate expres- 
sion of target proteins. This mechanism is in accord with 
the conceptual model of the present invention addressing 
the differentiation enigma, described hereinabove with 
specific reference to Figs. 6-7. 
[0263] Reference is now made to Fig. 18 which is a block diagram 
illustrating an overview of a methodology for finding novel 
oligonucleotides and operon-like polynucleotides of the 
present invention, and their respective functions. Accord- 
ing to a preferred embodiment of the present invention, 
the methodology to finding novel oligonucleotides of the 
present invention and their function comprises of the fol- 
lowing major steps: First, FIND GAM OLIGONUCLEOTIDES 
146 is used to detect , oligonucleotide of the novel group 



of oligonucleotide of the present invention, referred to 
here as GAM oligonucleotide. GAM oligonucleotides are 
located and their function elicited by detecting target pro- 
teins they bind and the function of those target proteins, 
as described hereinabove with reference to Figs. 9 - 15. 
Next, FIND GR POLYNUCLEOTIDES 147 is used to detect 
polynucleotide of a novel group of operon-like polynu- 
cleotide of the present invention, referred to here as GR 
polynucleotide. GR polynucleotides are located, by locat- 
ing clusters of proximally located GAM oligonucleotide, 
based on the previous step. Consequently, FIND HIERAR- 
CHY OF GR POLYNUCLEOTIDES 148 elicits the hierarchy of 
GR and GAM: binding sites for non-protein-binding GAM 
oligonucleotide comprised in each GR polynucleotide 
found are sought adjacent to other GR polynucleotides. 
When found, such a binding site indicates that the con- 
nection between the GAM and the GR the expression of 
which it modulates, and thus the hierarchy of the GR 
polynucleotides and the GAM oligonucleotides they com- 
prise. Lastly, DEDUCE FUNCTION OF HIGH GR POLYNU- 
CLEOTIDES AND GAM OLIGONUCLEOTIDES 149 is used to 
deduce the function of GR polynucleotides and GAM 
oligonucleotides which are ^high^ in the hierarchy, i.e. 



GAM oligonucleotides which modulate expression of other 
GR polynucleotides rather than directly modulating ex- 
pression of target proteins. A preferred approach is as 
follows: The function of protein-modulating GAM 
oligonucleotides is deducible from the proteins which they 
modulate, provided that the function of these target pro- 
teins is known. The function of ^ higher^ GAM oligonu- 
cleotides may be deduced by comparing the function of 
protein-modulating GAM oligonucleotides with the hierar- 
chical relationships by which the ^ higher^ GAM oligonu- 
cleotides are connected to the protein-modulating GAM 
oligonucleotides. For example, given a group of several 
protein-modulating GAM oligonucleotides which collec- 
tively cause a protein expression pattern typical of a cer- 
tain cell-type, then a "higher" GAM oligonucleotide is 
sought which modulates expression of GR polynucleotides 
which perhaps modulate expression of other GR polynu- 
cleotides which eventually modulate expression of the 
given group of protein- modulating GAM oligonucleotide. 
The "higher" GAM oligonucleotide found in this manner is 
taken to be responsible for differentiation of that cell- 
type, as per the conceptual model of the invention de- 
scribed hereinabove with reference to Fig. 6. 



[0264] Reference is now made to Fig. 19 wliicli is a blocl< diagram 
illustrating different utilities of oligonucleotide of the 
novel group of oligonucleotides of the present invention 
referred to here as CAM oligonucleotides and GR polynu- 
cleotides. The present invention discloses a first plurality 
of novel oligonucleotides referred to here as GAM 
oligonucleotides and a second plurality of operon-like 
polynucleotides referred to here as GR polynucleotides 
each of the GR polynucleotide encoding a plurality of GAM 
oligonucleotides The present invention further discloses a 
very large number of known target genes, which are 
bound by, and the expression of which is modulated by 
each of the novel oligonucleotides of the present inven- 
tion. Published scientific data referenced by the present 
invention provides specific, substantial, and credible evi- 
dence that the above mentioned target genes modulated 
by novel oligonucleotides of the present invention, are as- 
sociated with Alzheimers disease. Specific novel oligonu- 
cleotides of the present invention, target genes associated 
with Alzheimers diseases are described hereinbelow with 
reference to Tables 1 through 11. It is therefore appreci- 
ated that a function of GAM oligonucleotides and GR 
polynucleotides of the present invention is modulation of 



expression of target genes related to Alzheimers disease, 
and tliat therefore utilities of novel oligonucleotides of the 
present invention include diagnosis and treatment of 
Alzheimers disease. Fig. 19 describes various types of di- 
agnostic and therapeutic utilities of novel oligonucleotides 
of the present invention. A utility of novel oligonu- 
cleotides of the present invention is detection of GAM 
oligonucleotides and of GR polynucleotides. It is appreci- 
ated that since GAM oligonucleotides and polynucleotides 
modulate expression of disease related target genes, that 
detection of expression of GAM oligonucleotides in clini- 
cal scenarios associated with said diseases is a specific, 
substantial and credible utility. Diagnosis of novel 
oligonucleotides of the present invention may preferably 
be implemented by RNA expression detection techniques, 
including but not limited to biochips, as is well known in 
the art. Diagnosis of expression of oligonucleotides of the 
present invention may be useful for research purposes, in 
order to further understand the connection between the 
novel oligonucleotides of the present invention and the 
Alzheimers disease, for disease diagnosis and prevention 
purposes, and for monitoring disease progress. Another 
utility of novel oligonucleotides of the present invention is 



anti-GAM therapy, a mode of therapy which allows up 
regulation of the disease related target gene of a novel 
GAM oligonucleotide of the present invention, by lowering 
levels of the novel GAM oligonucleotide which naturally 
inhibits expression of that target gene. This mode of 
therapy is particularly useful with respect to target genes 
which have been shown to be under-expressed in associ- 
ation with Alzheimers disease. Anti-GAM therapy is fur- 
ther discussed hereinbelow with reference to Figs. 20A 
and 20B. A further utility of novel oligonucleotides of the 
present invention is GAM replacement therapy, a mode of 
therapy which achieves down regulation of Alzheimers re- 
lated target gene of a novel GAM oligonucleotide of the 
present invention, by raising levels of the GAM which nat- 
urally inhibits expression of that target gene. This mode 
of therapy is particularly useful with respect to target 
genes which have been shown to be over-expressed in 
Alzheimers disease. GAM replacement therapy involves 
introduction of supplementary GAM products into a cell, 
or stimulation of a cell to produce excess GAM products. 
GAM replacement therapy may preferably be achieved by 
transfecting cells with an artificial DNA molecule encoding 
a GAM which causes the cells to produce the GAM prod- 



uct, as is well known in the art. Yet a further utility of 
novel oligonucleotides of the present invention is modi- 
fied GAM therapy. Disease conditions are likely to exist, in 
which a mutation in a binding site of a GAM RNA prevents 
natural GAM RNA to effectively bind inhibit a disease re- 
lated target gene, causing up regulation of that target 
gene, and thereby contributing to the disease pathology. 
In such conditions, a modified GAM oligonucleotides is 
designed which effectively binds the mutated GAM bind- 
ing site, i.e. is an effective anti-sense of the mutated GAM 
binding site, and is introduced in disease effected cells. 
Modified GAM therapy is preferably achieved by transfect- 
ing cells with an artificial DNA molecule encoding the 
modified GAM which causes the cells to produce the mod- 
ified GAM product, as is well known in the art. An addi- 
tional utility of novel GAM of the present invention is in- 
duced cellular differentiation therapy. As aspect of the 
present invention is finding oligonucleotides which deter- 
mine cellular differentiation, as described hereinabove 
with reference to Fig. 18. Induced cellular differentiation 
therapy comprises transfection of cell with such GAM 
oligonucleotides thereby determining their differentiation 
as desired. It is appreciated that this approach may be 



widely applicable, inter alia as a means for auto trans- 
plantation harvesting cells of one cell-type from a patient, 
modifying their differentiation as desired, and then trans- 
planting them back into the patient. It is further appreci- 
ated that this approach may also be utilized to modify cell 
differentiation in vivo, by transfecting cells in a genetically 
diseased tissue with a cell-differentiation determining 
GAM thus stimulating these cells to differentiate appropri- 
ately. 

[0265] Reference is now made to Figs. 20A and 20B, simplified 
diagrams which when taken together illustrate anti-GAM 
therapy mentioned hereinabove with reference to Fig. 19. 
A utility of novel GAMs of the present invention is anti- 
GAM therapy, a mode of therapy which allows up regula- 
tion of Alzheimers disease related target gene of a novel 
GAM of the present invention, by lowering levels of the 
novel GAM which naturally inhibits expression of that tar- 
get gene. Fig. 20A shows a normal GAM inhibiting trans- 
lation of a target gene of GAM RNA by binding to a BIND- 
ING SITE found in an untranslated region of GAM TARGET 
RNA, as described hereinabove with reference to Fig. 8. 

[0266] Fig. 20B shows an example of anti-GAM therapy. ANTI- 
GAM RNA is short artificial RNA molecule the sequence of 



which is an anti-sense of CAM RNA. Anti-GAI\/I treatment 
comprises transfecting diseased cells with ANTI-GAM 
RNA, or with a DNA encoding thereof. The ANTI-CAM RNA 
binds the natural GAM RNA, thereby preventing binding of 
natural GAM RNA to its BINDING SITE. This prevents natu- 
ral translation inhibition of GAM TARGET RNA by GAM 
RNA, thereby up regulating expression of GAM TARGET 
PROTEIN. 

[0267] It is appreciated that anti-GAM therapy is particularly use- 
ful with respect to target genes which have been shown to 
be under-expressed in Alzheimers disease. Furthermore, 
anti-GAM therapy is particularly useful, since it may be 
used in situations in which technologies known in the art 
as RNAI and siRNA can not be utilized. As in known in the 
art, RNAI and siRNA are technologies which offer means 
for artificially inhibiting expression of a target protein, by 
artificially designed short RNA segments which bind com- 
plementarily to mRNA of said target protein. However, 
RNAi and siRNA can not be used to directly up regulate 
translation of target proteins. 

[0268] Reference is now made to Fig. 21A, which is a histogram 
representing the distribution of known mlRNA oligonu- 
cleotide s and that of hairpin structures extracted from 



expressed genome sequences with respect to their HAIR- 
PIN DETECTOR score. The known miRNA oligonucleotide s 
set is taken from REAM database, Release2.1 and include 
440 mlRNA oligonucleotide s from H.sapienas, 
M.musculus, C.elegans ,C.brigassae and D.melanogaster. 
Folding of expressed genome sequences taken from pub- 
lic databases of ESTs (Unigene-NCBI and TIGR) identified 
342,882 hairpin structures.~154,000 out of the 342,882 
hairpin structures did not pass the filter of being identi- 
fied as hairpins in several secondary structure folding ver- 
sions of the given genomic sequence , as described here- 
inabove with reference to Fig.l2B, and hence did not re- 
ceive a Hairpin detector score. Furthermore, ~133,000 
hairpin structures did not pass the filter of minimum 
score of the DICER-CUT LOCATION DETECTOR 116 (Fig. 9) 
(those ~287,000 hairpin structures are not represented in 
the histogram). Hairpin structures are considered as 
miRNA -like precursor oligonucleotides here referred to 
as GAM oligonucleotide, if their Hairpin detector score is 
above 0.3. Thus, the GAM oligonucleotides set is com- 
prised of ~40,000 hairpin structures, of those ~5100 re- 
ceived a high Hairpin detector score (>=0.7). These are 
much higher numbers than those of the known miRNA 



oligonucleotide s and of the upper bound of ~255 human 
miRNA oligonucleotide s, estimated by Bartel et al 
(Science,299,1540, March 2003). Of the reference set that 
pass the above filter (408/440), 284 (69%) received a high 
Hairpin detector score (>=0.7). 

[0269] Reference is now made to Fig. 2 IB, which is a table sum- 
marizing laboratory validation results that validate efficacy 
of the BIOINFORMATIC OLIGONUCLEOTIDE DETECTION 
ENGINE 100 ( Fig. 9). In order to assess efficacy of the 
BIOINFORMATIC OLIGONUCLEOTIDE DETECTION ENGINE 
100, novel oligonucleotides predicted thereby are prefer- 
ably divided into 4 detection accuracy groups (first col- 
umn), designated A through D, ranking GAMS from the 
most probable GAMs to the least probable GAMs, using 
the scores of HAIRPIN DETECTOR 114 (Fig. 9) and DICER- 
CUT LOCATION DETECTOR 116 (Fig. 9) as follows: 

[0270] Group A: The score of the HAIRPIN- DETECTOR is above 
0.7, the overall score of the two-phased predictor is 
above 0.55, and the score of the second phase of the 
two-phased predictor is above 0.75, or the score of the 
EDIT-DISTANCE predictor is equal or above 17. In this 
group, one Dicer cut location is predicted for each hairpin. 
Group B: The score of the HAIRPIN -DETECTOR is above 



0.5, the overall score of the two-phased predictor is 
above 0.55, and the hairpin is not in group A. Group C: 
The score of the HAIRPIN- DETECTOR is between 0.4 and 
0.5, and the overall score of the two-phased predictor is 
above 0.55. Group D: The score of the HAIRPIN-DETEC- 
TOR is between 0.3 and 0.4, and the overall score of the 
two-phased predictor is above 0.55. In groups B, C and D, 
if the score of the second phase of the two-phased pre- 
dictor is above 0.75, one Dicer cut location is predicted 
for each hairpin, otherwise both sides of the double 
stranded window are given as output, and are examined 
in the lab or used for binding site search. The groups are 
mutually exclusive, i.e. in groups A, C and D all hairpins 
score less than 17 in the EDIT-DISTANCE predictor. 
[0271] It is appreciated that the division into groups is not ex- 
haustive: 410 of the 440 published hairpins (second col- 
umn), and 1419 of the 1708 novel GAMs, belong to one of 
the groups. An indication of the real performance of the 
two-phased predictor in the presence of background hair- 
pins is given by the column 'precision on hairpin mixture' 
(third column). The precision on hairpin mixture is com- 
puted by mixing the published mlRNA hairpins with back- 
ground hairpins in a ratio of 1:4 and taking as a working 



assumption that they are hairpins not carrying a 'diced' 
miRNA.-lil<e oligonucleotide This is a strict assumption, 
since some of these background hairpins may indeed 
contain 'diced' miRNAs-like oligonucleotide, while in this 
column they are all counted as failures 

[0272] Sample novel bioinformatically predicted human GAMs of 
each of these groups are sent to the laboratory for valida- 
tion (fourth column), and the number (fifth column) and 
percent (sixth column) of successful validation of pre- 
dicted human GAM is noted for each of the groups, as 
well as overall (bottom line). The number of novel VAM 
genes explicitly specified by present invention belonging 
to each of the four groups is noted (seventh column). 

[0273] It is appreciated that the present invention comprises 

1419 novel GAM oligonucleotides, which fall into one of 
these four detection accuracy groups, and that the BIOIN- 
FORMATIC OLIGONUCLEOTIDE DETECTION ENGINE 100 ( 
Fig. 9) is substantiated by a group of 52 novel human 
GAM oligonucleotides validated by laboratory means, out 
of 168 human GAM oligonucleotides which were tested in 
the lab, resulting in validation of an overall 31% accuracy. 
The top group demonstrated 37% accuracy. Pictures of 
test-results of specific human GAM oligonucleotides in 



the abovementioned four groups, as well as the method- 
ology used for validating the expression of predicted 
oligonucleotides are elaborated hereinbelow with refer- 
ence to Fig. 22. 

[0274] It is further appreciated that failure to detect a predicted 
GAM oligonucleotide in the lab does not necessarily indi- 
cate a mistaken bioinformatic prediction. Rather, it may 
be due to technical sensitivity limitation of the lab test, or 
because the GAM oligonucleotide s not expressed in the 
tissue examined, or at the development phase tested. 

[0275] It is still further appreciated that in general these findings 
are in agreement with the expected bioinformatic accu- 
racy, as describe hereinabove with reference to Fig. 13B: 
assuming 80% accuracy of the HAIRPIN DETECTOR 114 
and 80% accuracy of the DICER-CUT LOCATION DETECTOR 
116 and 80% accuracy of the lab validation, this would re- 
sult in 50% overall accuracy of the GAM oligonucleotide 
validated in the lab. 

[0276] Reference is now made to Fig. 22A which is a picture of 
laboratory results validating the expression of 43 novel 
genes detected by the BIOINFORMATIC OLIGONUCLEOTIDE 
DETECTION ENGINE 100 (Fig. 9). 

[0277] Reference is now made to Fig. 22A and Fig.22B which are 



pictures and a summary table of laboratory results vali- 
dating the expression of 43 novel human GAM oligonu- 
cleotides detected by the BIOINFORMATIC OLIGONU- 
CLEOTIDE DETECTION ENGINEIOO. In each row in Fig.22A, 
pictures of several oligonucleotides validated by hy- 
bridization of Polymerase Chain Reaction (PCR)-product 
southern-blots, are provided, each corresponding to a 
specific GAM oligonucleotides, as elaborated hereinbelow. 
To test our validation method, we used a reference set of 
8 known human miRNA oligonucleotide s , as blind test to 
our laboratory. These PCR-product hybridization pictures 
are designated 1 through 8 for the reference set known 
miRNA oligonucleotide s; and 9 through 51 for predicted 
GAM oligonucleotides. 
[0278] In each PGR hybridization picture, 2 lanes are seen: the 

test lane, designated "+" and the control lane, designated 

For convenience of viewing the results, all PCR- 
product hybridization pictures of Fig. 22A have been 
shrunk x4 vertically. It is appreciated that for each of the 
tested GAM oligonucleotides a clear hybridization band 
appears in the test ("+") lane, but not in the control ("-") 
lane. 

[0279] Specifically, Fig. 22A shows pictures of PCR-product hy- 



bridization validation by soutliern-blot, tlie methodology 
of wliicli is described liereinbelow, to the following novel 
human GAM oligonucleotides (RosettaGenomics Ltd. 
Nomenclature, 'A' and 'B' referred to the Dicer Cut Loca- 
tion as described hereinabelow with reference to the de- 
scription of large tables: 
[0280] (i)hsa-MIR-21;(2)hsa-MIR-27b; (3)hsa-MIR-186; 
(4)hsa-MIR-93; (5)hsa-MIR-26a ; (6)hsa-MIR-191; 
(7)hsa-MIR-31; (8)hsa-MIR-92; (9) GAM3418-A (later 
published by other researchers as hsa-MIR23); (10) 
GAM4426-A; (11) GAM281-A; (12) GAM7553-A; (13) 
GAM5385-A; (14) GAM2608-A; (15) GAM1032-A; (16) 
GAM3431-A; (17) GAM7933-A; (18) GAM3298-A.; (19) 
GAM7080-A; (20)GAM895-A.; (21) GAM3770.1; (22) 
GAM337162-A; (23) GAM 8678-A; (24) GAM2033-A; 
(25)GAM7776-A; (26) GAM8145-A; (27) GAM25-A; (28) 
GAM7352.1; (29) GAM337624-A; (30) GAM1479-A; (31) 
GAM2270-A; (32) GAM7591-A; (33) GAM8285-A; (34) 
GAM6773-A; (35) GAM336818-A; (36) GAM336487-A; 
(37) GAM337620-A; (38) GAM336809-A; (39) 
GAM5346-A; (40) GAM8554-A; (41)GAM2071-A; (42) 
GAM7957-A; (43) GAM391-A; (44) GAM6633-A; (45) 
GAM19; (46) GAM8358-A; (47) GAM3229-A; an)GAM 



7052-A; (49) GAM3027-A; (50) GAM21and (51) GAM 
oligonucleotide similar to mmu-MIR-30e. 
[0281] The next validated GAM oligonucleotides are highly simi- 
lar or highly identical to known mouse-mlRNA oligonu- 
cleotides: GAM3027-A, similar to mmu-MIR-29c; GAM21, 
similar to mmu-MIR-130b; and GAM oligonucleotide 
which is highly similar to mmu-MIR-30e (picture number 
51). In addition to the PGR- product hybridization detec- 
tion, the following GAMs were cloned and sequenced : 
GAM3418-A , GAM5385-A, GAM1032-A, GAM3298-A, 
GAM7080-A, GAM1338-A, GAM7776-A, GAM25-A, 
GAM337624-A, GAM1479-A, GAM6773-A, 
GAM336818-A, GAM336487-A, GAM337620-A, 
GAM336809-A, GAM3027-A, GAM21, and GAM oligonu- 
cleotide similar to mmu-MIR-30e (picture number 51). 
Furthermore, the following GAM oligonucleotides were se- 
quenced directly from the ligation reaction by the method 
described hereinbelow under LIGATION-PCR DIAGNOSTIC 
METHOD: GAM4426-A, GAM7553-A, GAM2270-A, and 
GAM7591-A. 

[0282] In order to validate the expression of predicted novel GAM 
and assuming that these novel GAM oligonucleotides are 
probably expressed at low concentrations, a PGR product 



cloning approach was set up through the following strat- 
egy: two types of cDNA libraries designated "One tailed" 
and "Ligation" were prepared from frozen HeLa SlOO ex- 
tract (4c Biotech, Belgium) size fractionated RNA. Essen- 
tially, Total SlOO RNA was prepared through an SDSPro- 
teinase K incubation followed by an acid Phenol- 
Chloroform purification and Isopropanol precipitation. Al- 
ternatively, total HeLa RNA was also used as starting ma- 
terial for these libraries. 
[0283] Fractionation was done by loading up to 500g per YMIOO 
Amicon Microcon column (Millipore) followed by a 500g 
centrifugation for 40 minutes at 4C. Flow through 
"YM100"RNA consisting of about of the total RNA was 
used for library preparation or fractionated further by 
loading onto a YM30 Amicon Microcon column (Millipore) 
followed by a 13,500g centrifugation for 25 minutes at 
4C. Flowthrough "YM30" was used for library preparation 
as is and consists of less than 0.5% of total RNA. For the 
both the "ligation" and the "One-tailed" libraries, RNA was 
dephosphorilated and ligated to an RNA (lowercase)-DNA 
(UPPERCASE) hybrid 5"-phosphorilated, 3" idT blocked 
3"-adapter (5"-P-uuuAACCGCATCCTTCTC-idT-3" Dhar- 
macon # P-002045-01-05) (as elaborated in Elbashir et 



al., Genes Dev. 15:188-200 (2001)) resulting in ligation 
only of RNase III type cleavage products. 3"-Ligated RNA 
was excised and purified from a half 6%, half 13% poly- 
acrylamide gel to remove excess adapter with a Nanosep 
0.2M centrifugal device (Pall) according to instructions, 
and precipitated with glycogen and 3 volumes of Ethanol. 
Pellet was resuspended in a minimal volume of water. 
[0284] For the "ligation" library a DNA (UPPERCASE)- RNA 
(lowercase) hybrid 5"-adapter 
(5"-TACTAATACGACTCACTaaa-3" Dharmacon # P- 
002046-01-05) was ligated to the 3"-adapted RNA, re- 
verse transcribed with "EcoRI-RT": 
(5"-GACTAGCTGGAATTCAAGGATGCGGTTAAA-3"), PGR 
amplified with two external primers essentially as in El- 
bashir et al 2001 except that primers were "EcoRI-RT'and 
"PstI 

Fwd"(5"-CAGCCAACGCTGCAGATACGACTCACTAAA-3"). 
This PGR product was used as a template for a second 
round of PGR with one hemispecific and one external 
primer or with two hemispecific primers. 
[0285] For the "One tailed" library the 3"-Adapted RNA was an- 
nealed to 20pmol primer "EcoRI RT" by heating to 70G and 
cooling O.lG/sec to 30G and then reverse transcribed with 



Superscript II RT (According to instructions, Invitrogen) in 
a 201 volume for 10 alternating 5 minute cycles of 37C 
and 45C. Subsequently, RNA was digested with II 2M 
NaOH, 2mM EDTA at 65C for 10 minutes. cDNA was 
loaded on a polyacrylamide gel, excised and gel-purified 
from excess primer as above (invisible, judged by primer 
run alongside) and resuspended in 131 of water. Purified 
cDNA was then oligo-dC tailed with 400U of recombinant 
terminal transferase (Roche molecular biochemicals), II 
lOOM dCTP, II 15mM CoCI2, and 41 reaction buffer, to a 
final volume of 201 for 15 minutes at 37C. Reaction was 
stopped with 21 0.2M EDTA and 151 3M NaOAc pH 5.2. 
Volume was adjusted to 1501 with water. Phenol : Bro- 
mochloropropane 10:1 extracted and subsequently pre- 
cipitated with glycogen and 3 volumes of Ethanol. C-tailed 
cDNA was used as a template for PGR with the external 
primers 

"T3-PstBsg(G/l)18"(5"-AATTAACCCTCACTAAAGGCTGCAG 
GTGCAGGIGGGIIGGGIIGGGIIGN-3" where I stands for Ino- 
sine and N for any of the 4 possible deoxynucleotides), 
and with "EcoRI 

Nested"(5"-GGAATTGAAGGATGGGGTTA-3"). This PGR 
product was used as a template for a second round of PGR 



with one hemispecific and one external primer or witli two 
hemispecific primers. 
[0286] Hemispecific primers were constructed for each predicted 
CAM RNA oligonucleotide by an in-house program de- 
signed to choose about half of the 5" or 3" sequence of 
the GAM RNA corresponding to a TM of about 30-34C 
constrained by an optimized 3" clamp, appended to the 
cloning adapter sequence (for "One-tailed" libraries 
5"-GGNNGGGNNG on the 5" end of the GAM RNA , or TT- 
TAACCGCATC-3" on the 3" end of the GAM RNA. For "Liga- 
tion" libraries the same 3" adapter and 
5"-CGACTCACTAAA on the 5" end). Consequently, a fully 
complementary primer of a TM higher than 60C was cre- 
ated covering only one half of the GAM RNA sequence 
permitting the unbiased elucidation by sequencing of the 
other half. 

[0287] CONFIRMATION OF GAM OLIGONUCLEOTIDE SEQUENCE 
AUTHENTICITY OF PGR PRODUCTS: 

[0288] SOUTHERN BLOT: PGR- product sequences were confirmed 
by southern blot (Southern E.M., Biotechnology, 
1992,24:122-39 (1975)) and hybridization with DNA 
oligonucleotide probes synthesized against predicted 
GAM RNAs oligonucleotides. Gels were transferred onto a 



Biodyne PLUS 0.45m, (Pall) positively charged nylon mem- 
brane and UV cross-linked. Hybridization was performed 
overnight with DIG-labeled probes at 420C in DIG Easy- 
Hyb buffer (Roche). Membranes were washed twice with 
2xSSC and 0.1% SDS for 10 min. at 420C and then washed 
twice with O.SxSSC and 0.1% SDS for 5 min at 420C. The 
membrane was then developed by using a DIG lumines- 
cent detection kit (Roche) using anti-DIG and CSPD reac- 
tion, according to the manufacturer's protocol. All probes 
were prepared according to the manufacturers (Roche 
Molecular Biochemicals) protocols: Digoxigenin (DIG) la- 
beled antisense transcripts was prepared from purified 
PGR products using a DIG RNA labeling kit with 73 RNA 
polymerase. DIG labeled PGR was prepared by using a DIG 
PGR labeling kit. 3"-DIG-tailed oligo ssDNA antisense 
probes, containing DIG-dUTP and dATP at an average tail 
length of 50 nucleotides were prepared from lOOpmole 
oligonucleotides with the DIG Oligonucleotide Labeling 
Kit. 

[0289] GLONE-SEQUENGING: PGR products were inserted into 
pGEM-T (Promega) or pTZ57 (MB! Fermentas), trans- 
formed into competent JM109 E. coll (Promega) and sown 
on LB-Amp plates with IPTG/Xgal. White and light-blue 



colonies were transferred to duplicate gridded plates, one 
of which was blotted onto a membrane (Biodyne Plus, Pall) 
for hybridization with DIG tailed oligo probes (according 
to instructions, Roche) corresponding to the expected 
GAM. Plasmid DNA from positive colonies was sequenced. 
[0290] LIGATION-PCR DIAGNOSTIC METHOD: To further validate 
predicted GAM PGR product sequence derived from hemi- 
primers, a PGR based diagnostic technique was devised to 
amplify only those products containing also at least two 
additional nucleotides of the non hemi-primer defined 
part of the predicted GAM RNA oligonucleotide. In 
essence, a diagnostic primer was designed so that its 3" 
end, which is the specificity determining side, was identi- 
cal to the desired GAMRNA oligonucleotide, 2-10 nu- 
cleotides (typically 4-7, chosen for maximum specificity) 
further into its 3" end than the nucleotide stretch primed 
by the hemi-primer. The hemi-primer PGR product was 
first ligated into a T-cloning vector (pTZ57/T or pGEM-T) 
as described herinabove. The ligation reaction mixture 
was used as template for the diagnostic PGR under strict 
annealing conditions with the new diagnostic primer in 
conjunction with a general plasmid-homologous primer, 
resulting in a distinct ~200 base-pair product. This PGR 



product can be directly sequenced, permitting tlie eluci- 
dation of the remaining nucleotides up to the 3" of the 
mature GAM RNA oligonucleotide adjacent to the 3" 
adapter. Alternatively, following analysis of the diagnostic 
PGR reaction on an agarose gel, positive ligation reactions 
(containing a band of the expected size) were transformed 
into E. coll. Using this same diagnostic technique and as 
an alternative to screening by Southern-blot colony- 
hybridization, transformed bacterial colonies were 
screened by colony-PGR (Gussow,D. and Glackson,T, Nu- 
cleic Acids Res. 17: 4000 (1989)) prior to plasmid purifi- 
cation and sequencing. 
[0291] Reference is now made to Fig. 22B which is a table sum- 
marizing laboratory results which validate the expression 
of 8 known human mlRNA oligonucleotide s and 43 novel 
GAM oligonucleotides detected by the BIOINFORMATIG 
OLIGONUGLEOTIDE DETEGTION ENGINE 100. The table 
gives additional information on the above GAM oligonu- 
cleotides and contains the following fields:: NUMBER: refer 
to the hybridization picture number of Fig.22A; NAME: in- 
dicate the known MIR name for the reference set or GAM's 
name as given by RosettaGenomices nomenclature 
method; SEQUENGE: 5' to 3' sequence of the mature , 



'diced' oligonucleotide; SEQUENCED: '+' indicates a valida- 
tion of the GAM RNA sequence by sequencing procedure 
as described hereinabove with reference to Fig.22A. 
[0292] Reference is now made to Fig. 23A, which is a schematic 
representation of a novel human GR polynucleotide herein 
designated GR12731 (RosettaGenomics Ltd. nomencla- 
ture), located on chromosome 9, comprising 2 known hu- 
man MIR genes - MIR24 and MIR23, and 2 novel GAM 
oligonucleotides, herein designated GAM22 and GAM116, 
all marked by solid black boxes. Fig. 23A also schemati- 
cally illustrates 6 non-GAM hairpin sequences, and one 
non-hairpin sequence, all marked by white boxes, and 
serving as negative controls. By "non-GAM hairpin se- 
quences" is meant sequences of a similar length to known 
MIR PRECURSOR sequences, which form hairpin secondary 
folding pattern similar to MIR PRECURSOR hairpins, and 
yet which are assessed by the BIOINFORMATIC OLIGONU- 
CLEOTIDE DETECTION ENGINE 100 not to be valid GAM 
PRECURSOR hairpins. It is appreciated that Fig. 23A is a 
simplified schematic representation, reflecting only the 
order in which the segments of interest appear relative to 
one another, and not a proportional distance between the 
segments. 



[0293] Reference is now made to Fig. 23B, wliicli is a scliematic 
representation of secondary folding of each of the M\Rs 
and CAMS of CR GR12731 MIR24, IV1IR23, GAM22 and 
CAIS/1116, and of the negative control non-CAM hairpins, 
herein designated N2, N3, N116, N4, N6 and N7. NO is a 
non-hairpin control, of a similar length to that of known 
MIR PRECURSOR hairpins. It is appreciated that the nega- 
tive controls are situated adjacent to and in between real 
MIR genes and GAM predicted oligonucleotide and 
demonstrates similar secondary folding patterns to that of 
known MIRs and GAMs. 

[0294] Reference is now made to Fig. 23C, which is a picture of 
laboratory results of a PGR test upon a YMIOO "liga- 
tion"-library, utilizing specific primer sets directly inside 
the boundaries of the hairpins. Due to the nature of the li- 
brary the only PGR amplifiable products can result from 
RNaselll type enzyme cleaved RNA, as expected for legiti- 
mate hairpin precursors presumed to be produced by 
DROSHA (Lee et al. Nature 425 415-419, 2003). Fig 23G 
demonstrates expression of hairpin precursors of known 
MIR genes - MIRhsa-23 and MIRhsa-24, and of novel 
bioinformatically detected GAM22 and GAM 116 hairpins 
predicted bioinformatically by a system constructed and 



operative in accordance with a preferred embodiment of 
the present invention. Fig. 23C also shows that none of 
the 7 controls (6 hairpins designated N2, N3, N23, N4, N6 
and N7 and 1 non-hairpin sequence designated NO) were 
expressed. N 116 is a negative control sequence partially 
overlapping GAM116. 
[0295] In the picture, test lanes including template are desig- 
nated "+" and the control lane is designated It is ap- 
preciated that for each of the tested hairpins, a clear PGR 
band appears in the test ("+") lane, but not in the control 
("-") lane. 

[0296] Figs. 23A through 23C, when taken together validate the 
efficacy of the bioinformatic oligonucleotide detection en- 
gine in: (a) detecting known MIR genes; (b) detecting novel 
CAM PRECURSOR hairpins which are found adjacent to 
these MIR genes, and which despite exhaustive prior bio- 
logical efforts and bioinformatic detection efforts, went 
undetected; (c) discerning between CAM (or MIR) PRECUR- 
SOR hairpins, and non-CAM hairpins. 

[0297] It is appreciated that the ability to discern CAM-hairpins 
from non-CAM-hairpins is very significant in detecting 
CAM oligonucleotide since hairpins in general are highly 
abundant in the genome. Other MIR prediction programs 



have not been able to address this challenge successfully. 

[0298] Reference is now made to Fig. 24A which is an annotated 
sequence of an EST comprising a novel GAM oligonu- 
cleotides detected by the oligonucleotide detection sys- 
tem of the present invention. Fig. 24A shows the nu- 
cleotide sequence of a known human non-protein coding 
EST (Expressed Sequence Tag), identified as EST72223. 
The EST72223 clone obtained from TIGR database 
(Kirkness and Kerlavage, 1997) was sequenced to yield the 
above 705bp transcript with a polyadenyl tail. It is appre- 
ciated that the sequence of this EST comprises sequences 
of one known mlRNA oligonucleotide , identified as hsa- 
MIR98, and of one novel GAM oligonucleotide referred to 
here as GAM25, detected by the BIOINFORMATIC 
OLIGONUCLEOTIDE DETECTION ENGINE 100 (Fig. 9) of the 
present invention. 

[0299] The sequences of the precursors of the known MIR98 and 
of the predicted GAM25 are precursor in bold, the se- 
quences of the established mlRNA 98 and of the predicted 
miRNA-like oligonucleotide GAM25 are underlined. 

[0300] Reference is now made to Figs. 24B, 24C and 24D that are 
pictures of laboratory results, which when taken together 
demonstrate laboratory confirmation of expression of the 



bioinformatically detected novel oligonucleotide of Fig. 
24A. In two parallel experiments, an enzymatically syn- 
thesized capped, EST72223 RNA transcript, was incubated 
with Hela SlOO lysate for 0 minutes, 4 hours and 24 
hours. RNA was subsequently harvested, run on a dena- 
turing poly aery I amide gel, and reacted with a 102nt and a 
145nt antisense MIR98 and GAM25 precursor transcript 
probes respectively. The Northern blot results of these 
experiments demonstrated processing of EST72223 RNA 
by Hela lysate (lanes 2-4, in 24B and 24C), into ~80bp 
and ~22bp segments, which reacted with the MIR98 pre- 
cursor probe (24B), and into ~100bp and ~24bp seg- 
ments, which reacted with the GAM25 precursor probe 
(24C). These results demonstrate the processing of 
EST72223 by Hela lysate into MIR98 precursor and GAM25 
precursor. It is also appreciated from Fig. 24C (lane 1) that 
Hela lysate itself reacted with the GAM25 precursor probe, 
in a number of bands, including a ~100bp band, indicat- 
ing that GAM25-precursor is endogenously expressed in 
Hela cells. The presence of additional bands, higher than 
lOObp in lanes 5-9 probably corresponds to the presence 
of nucleotide sequences in Hela lysate, which contain the 
GAM25 sequence. 



[0301] In addition, in order to demonstrate tlie l<inetics and 

specificity of tlie processing of M\R9S and GAIVIZS precur- 
sors into tlieir respective mature, 'diced' segments, tran- 
scripts of l\/IIR98 and of tlie bioinformatically predicted 
CAI\/I25 precursors were similarly incubated with Hela 
SlOO lysate, for 0 minutes, 30 minutes, 1 hour and 24 
hours, and for 24 hours with the addition of EDTA, added 
to inhibit Dicer activity, following which RNA was har- 
vested, run on a polyacrylamide gel and reacted with 
MIR98 and GAM25 precursor probes. Capped transcripts 
were prepared for in-vitro RNA cleavage assays with T7 
RNA polymerase including a m7G(5')ppp(5')G-capping re- 
action the Message Machine kit (Ambion). Purified PGR 
products were used as template for the reaction. These 
were amplified for each assay with specific primers con- 
taining a T7 promoter at the 5" end and a T3 RNA poly- 
merase promoter at the 3"end. Gapped RNA transcripts 
were incubated at 30G in supplemented, dialysis concen- 
trated, Hela SlOO cytoplasmic extract (4G Biotech, Seneffe, 
Belgium). The Hela SlOO was supplemented by dialysis to 
a final concentration of 20mM Hepes, lOOmM KGI, 2.5mM 
MgGI2 , 0.5mM DTT, 20% glycerol and protease inhibitor 
cocktail tablets (Complete mini Roche Molecular Biochem- 



icals). After addition of all components, final concentra- 
tions were lOOmM capped target RNA, 2mM ATP, 0.2mM 
GTP, 500U/ml RNasin, 25g/ml creatine kinase, 25mM cre- 
atine phosphate, 2.5mM DTT and 50%S100 extract. Pro- 
teinase K, used to enhance Dicer activity (Zhang et al., 
EMBOJ. 21, 5875-5885 (2002)) was dissolved in 50mM 
Tris-HCI pH 8, 5mM CaCI2, and 50% glycerol, was added 
to a final concentration of 0.6 mg/ml. Cleavage reactions 
were stopped by the addition of 8 volumes of proteinase K 
buffer (200Mm Tris-Hcl, pH 7.5, 25m M EDTA, 300mM 
NaCI, and 2% SDS) and incubated at 65C for 15min at dif- 
ferent time points (0, 0.5, 1, 4, 24h) and subjected to 
phenol/chloroform extraction. Pellets were dissolved in 
water and kept frozen. Samples were analyzed on a seg- 
mented half 6%, half 13% polyacrylamide 1XTBE-7M Urea 
gel. 

[0302] The Northern blot results of these experiments demon- 
strated an accumulation of a ~22bp segment which re- 
acted with the MIR98 precursor probe, and of a ~24bp 
segment which reacted with the GAM25 precursor probe, 
over time (lanes 5-8). Absence of these segments when 
incubated with EDTA (lane 9), which is known to inhibit 
Dicer enzyme (Zhang et al., 2002), supports the notion 



that the processing of MIR98 and GAM25 precursors into 
their ' diced' segments is mediated by Dicer enzyme, 
found in Hela lysate. The molecular sizes of EST72223, 
MIR-98 and GAM25 and their corresponding precursors 
are indicated by arrows. 

[0303] Fig. 24D present Northern blot results of same above ex- 
periments with GAM25 probe (24nt). The results clearly 
demonstrated the accumulation of mature GAM25 
oligonucleotide after 24 h. 

[0304] To validate the identity of the band shown by the lower 
arrow in figs. 24G and 24D, a RNA band parallel to a 
marker of 24 base was excised from the gel and cloned as 
in Elbashir et al (2001) and sequenced. 90 clones corre- 
sponded to the sequence of mature GAM25 oligonu- 
cleotide .three corresponded to GAM25* (the opposite 
arm of the hairpin with a 1-3 nucleotide 3" overhang) and 
two to the hairpin-loop. 

[0305] GAM25 was also validated endogenously by sequencing 
from both sides from a HeLa YMIOO total-RNA "ligation" 
libraries, utilizing hemispecific primers as described in 
Fig. 22. 

[0306] Taken together, these results validate the presence and 
processing of a novel MIR-like oligonucleotide, GAM25, 



which was predicted bioinformatically. The processing of 
this novel CAM oligonucleotide product, by Hela lysate 
from EST72223, through its precursor, to its final form 
was similar to that observed for known mlRNA oligonu- 
cleotide, MIR98. 

[0307] Transcript products were 705nt (EST72223), 102nt (MIR98 
precursor), 125nt (GAM25 precursor) long. EST72223 was 
PGR amplified with T7-EST 72223 forward primer: 
5"-TAATACGACTCACTATAGGCCCTTATTAGAGGATTCTGC 
T-3" and T3-EST72223 reverse 

primer:"-AATTAACCCTCACTAAAGGI II II III ICCTGAGA 
CAGAGT-3".MIR98 was PGR amplified using EST72223 as 
a template with T7MIR98 forward primer: 
5-"TAATACGACTCACTATAGGGTGAGGTAGTAAGTTGTATT 
GTT-3"and T3MIR98 reverse primer: 
5 "- AATTAACCCTCACTAAAGGG AAAGTAGTAAGTTGTATAG 
TT-3".GAM25 was PGR amplified using EST72223 as a 
template with GAM25 forward primer: 
5"-GAGGGAGGAGAATTGGTTGA-3" and T3-EST72223 re- 
verse 

primer:5"-AATTAAGGGTGAGTAAAGGGGTGAGACAGAGTCT 
TGGTG-3". 

[0308] It is appreciated that the data presented in Figs. 24A, 24B, 



24C and 24D when taken together validate the function of 
the bioinformatic oligonucleotide detection engine 100 of 
Fig. 9. Fig. 24A shows a novel GAM oligonucleotide bioin- 
formatically detected by the BIOINFORMATIC OLIGONU- 
CLEOTIDE DETECTION ENGINE 100, and Figs. 24C and 24D 
show laboratory confirmation of the expression of this 
novel oligonucleotide. This is in accord with the engine 
training and validation methodology described herein- 
above with reference to Fig. 10. 
DETAILED DESCRIPTION OF LARGE TABLES 

[0309] Table 1 comprises data relating the SEQ ID NO of GAM 
RNA oligonucleotides of the present invention to their 
corresponding GAM NAME, and contains the following 
fields: GAM SEQ-ID: GAM SEQ ID NO, as in the Sequence 
Listing; GAM NAME: Rosetta Genomics Ltd. nomenclature 
(see below); GAM RNA SEQUENCE: Sequence (5^ to 3^) of 
the mature, "diced" GAM RNA; GAM POS: Dicer cut loca- 
tion (see below); and 

[0310] Table 2 comprises detailed textual description according 
to the description of Fig. 8 of each of a plurality of novel 
GAM oligonucleotide of the present invention, and con- 
tains the following fields: GAM NAME: Rosetta Genomics 
Ltd. nomenclature (see below);PRECUR SEQ-ID:GAM pre- 



cursor Seq-ID, as in the Sequence Listing; PRECURSOR SE- 
QUENCE: Sequence (5^ to 3^) of tlie CAM precursor ; GAM 
DESCRIPTION: Detailed description of GAM oligonucleotide 
with reference to Fig. 8; and 
[0311] Table 3 comprises data relating to the source and location 
of novel GAM oligonucleotides of the present invention, 
and contains the following fields: GAM NAME: Rosetta Ge- 
nomics Ltd. nomenclature (see below); PRECUR SEQ-ID: 
GAM precursor SEQ ID NO, as in the Sequence Listing; OR- 
GANISM: Abbreviated (hsa = Homo sapiens); CHR: Chro- 
mosome encoding the GAM oligonucleotide; STRAND: Ori- 
entation on the chromosome, '+' for the plus strand,'-' for 
the minus strand; CHR- START OFFSET Start offset of GAM 
precursor sequence on the chromosome; CHR-END OFF- 
SET: End offset of GAM precursor sequence on the chro- 
mosome; SOURCE_REF-ID: Accession number of source 
sequence; and 

[0312] Table 4 comprises data relating to GAM precursors of 

novel GAM oligonucleotides of the present invention, and 
contains the following fields: GAM NAME: Rosetta Ge- 
nomics Ltd. nomenclature (see below); PRECUR SEQ-ID: 
GAM precursor Seq-ID, as in the Sequence Listing; PRE- 
CURSOR-SEQUENCE: Sequence (5^ to 3^) of the GAM pre- 



cursor; GAM FOLDED PRECURSOR RNA: Schematic repre- 
sentation of tlie CAM folded precursor, beginning 5" end 
(beginning of upper row) to 3^ end (beginning of lower 
row), where the hairpin loop is positioned at the right part 
of the draw; and 

[0313] Table 5 comprises data relating to GAM oligonucleotides 
of the present invention, and contains the following fields: 
GAM NAME: Rosetta Genomics Ltd. nomenclature (see be- 
low); GAM RNA SEQUENCE: Sequence (5^ to 3^) of the ma- 
ture, ^diced^ GAM RNA ; PRECUR SEQ-ID : GAM precursor 
Seq-ID, as in the Sequence Listing; SOURCE_REF-ID: ac- 
cession number of the source sequence; GAM POS: Dicer 
cut location (see below); and 

[0314] Table 6 comprises data relating SEQ ID NO of the GAM 
target gene binding site sequence to TARGET gene name 
and target binding site sequence, and contains the follow- 
ing fields: TARGET BINDING SITE SEQ-ID: Target binding 
site SEQ ID NO, as in the Sequence Listing; TARGET: GAM 
target gene name; TARGET BINDING SITE SEQUENCE: Nu- 
cleotide sequence (5^ to 3^) of the target binding site; 
and 

[0315] Table 7 comprises data relating to target genes and bind- 
ing sites of GAM oligonucleotides of the present inven- 



tion, and contains the following fields: GAM NAME: 
Rosetta Genomics Ltd. nomenclature (see below); GAM 
RNA SEQUENCE: Sequence (5' to 3^) of the mature, 
"diced" GAM RNA; TARGET: GAM target gene name; TAR- 
GET REF-ID: Target accession number (GenBank); UTR: 
Untranslated region of binding site/s (3" or 5"); TARGET 
BS-SEQ: Nucleotide sequence (5' to 3^) of the target 
binding site; BINDING-SITE-DRAW: Schematic representa- 
tion of the binding site, upper row represent 5" to 3" se- 
quence of the GAM RNA, lower row represent 3" to 5" se- 
quence of the target binding site; GAM POS: Dicer cut lo- 
cation (see below); and 

[0316] Table 8 comprises data relating to functions and utilities 
of novel GAM oligonucleotides of the present invention, 
and contains the following fields: GAM NAME: Rosetta Ge- 
nomics Ltd. nomenclature (see below);TARGET: GAM tar- 
get gene name; GAM RNA SEQUENCE: Sequence (5" to 3") 
of the mature, "diced" GAM RNA; GAM FUNCTION: De- 
scription of the GAM functions and utilities; GAM POS: 
Dicer cut location (see below);TAR DIS: Target Disease Re- 
lation Group (see below); and 

[0317] Table 9 comprises data of GAM target gene function ref- 
erences- Bibliography and contains the following fields: 



GAM NAME: Rosetta Genomics Ltd. nomenclature (see be- 
low); GAM RNA SEQUENCE: Sequence (5^ to 3^) of the ma- 
ture, ^ diced ^ GAM RNA; TARGET: GAM target gene name; 
REFERENCES: list of references related to the GAM target 
gene; GAM POS: Dicer cut location (see below); and 
[0318] Table 10 comprises data relating to novel GR (Genomic 
Record) polynucleotides of the present invention, and 
contains the following fields: GR NAME: Rosetta Genomics 
Ltd. nomenclature (see below); GR DESCRIPTION: Detailed 
description of a GR polynucleotide cluster, with reference 
to Fig. 16; and 

[0319] Table 11 comprises data relating to Alzheimers disease 
that GAM oligonucleotides are predicted to regulate the 
disease-associated genes. Each row is referred to a spe- 
cific disease, and list the GAM target genes related to the 
disease. The first row is a summary of ALL target genes 
associated in Alzheimer disease containing in the present 
invention. The second row is a subset of the first row and 
contains all GAM target genes found to bind to at least 
one validated GAM oligonucleotide. The table contains the 
following fields: ROW#: index of the row number; DISEASE 
NAME: name of the disease; TARGET GENES ASSOCIATED 
WITH ALZHEIMER: list of GAM target genes that are asso- 



ciated with the specified disease; and 
[0320] The following conventions and abbreviations are used in 
the tables: The nucleotide 'U' is represented as T in the 
tables, and 

[0321] GAM NAME or GR NAME are names for nucleotide se- 
quences of the present invention given by RosettaGe- 
nomics Ltd. nomenclature method. All GAMs/GRs are des- 
ignated by GAMx/GRx where x is a unique ID. 

[0322] SOURCE REF-ID: The accession number of expressed se- 
quences on which novel oligonucleotides were detected. 
The sequences are taken from the following published 
databases: (1) TIGR- "Tentative Human Consensus" (THC) 
(2) EST database-UNIGENE, NCBI. 

[0323] GAM POS is a position of the GAM RNA on the GAM PRE- 
CURSOR RNA sequence. This position is the Dicer cut lo- 
cation, 'A' indicates a probable Dicer cut location, 'B' indi- 
cates an alternative Dicer cut location. 

[0324] TAR DIS (Target Disease Relation Group) "A" indicates if the 
target gene is known to have a specific causative relation 
to Alzheimers disease, based on the OMIM database 
(Hamosh et al, 2002). It is appreciated that this is a partial 
classification emphasizing genes which are associated 
with "single gene" diseases etc. All GAM oligonucleotides 



of the present invention ARE associated witli Alzlieimers 
disease, although not all are necessary in 'A' status. 
[0325] All genomic sequences of the present invention as well as 
their chromosomal location and strand orientation are de- 
rived from sequences records of NCBI, Build33 database 
(April, 2003). 

[0326] It is appreciated by persons skilled in the art that the 

present invention is not limited by what has been particu- 
larly shown and described hereinabove. Rather the scope 
of the present invention includes both combinations and 
subcombinations of the various features described here- 
inabove as well as variations and modifications which 
would occur to persons skilled in the art upon reading the 
specifications and which are not in the prior art. 



